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An Automatic Recording Apparatus for the Study of Flow 
and Recovery in Metals 


By E. N. pa-C. ANDRADE anp A. J. KENNEDY 
University College, London 


MS. received 7th July 1949 


ABSTRACT. An apparatus is described which records continuously on photographic 
paper the extension-against-time curve of a metal wire creeping under stress. The length 
and time scales are recorded on the paper at the same time as the extension, so that deforma- 
tion of the paper during development and drying does not affect the accuracy of the record, 
which allows an extension of up to 10 cm. to be read to within 0:02 mm. The apparatus also 
automatically removes and restores, repeatedly if desired, the load at times which can be set 
before the start of the experiment. 


§1. INTRODUCTION 


s part of the work on the flow of metals which is being carried out in this 
laboratory it has been decided to investigate the effect of periods of 
recovery, during which the metal is not exposed to external stress, these 

periods being of varied length and interpolated at various stages of the flow. In 
such experiments it is very convenient to be able to record automatically the 
progressive extension of the metal against time, both when the stress is continu- 
ously applied and when it is removed for specified periods and then reapplied. 
The apparatus here described was designed for such automatic recording and 
control and has now been sufficiently used to establish that it answers its purpose 
well, 
§2. GENERAL DESCRIPTION 

The record is made on photographic paper wrapped round a drum which is 
driven at uniform speed by a phonic motor, suitably geared down. ‘The same 
motor drives another shaft which carries a timing disc furnished with contact 
studs. This shaft and the drum shaft are provided with maynetic clutches. 
The phonic motor also actuates a cam and switch which flashes a lamp, at short 
regular intervals, to provide a time scale on the paper. Horizontal regularly 
spaced lines, providing ascale of length, are laid down by the projection ofa vertical 
scale illuminated by a steady light source. 

The wire is hung vertically and the load applied through the :ntermediacy of 
a vertical rod attached to its lower end. ‘The wire can be relieved of stress by a 
device which grasps this rod and lifts it slightly, this device being actuated by the 
timing disc. At the sameinstant the lights which register the time and lengthscales 
and the instantaneous length of the wire are switched off and the drum 1s stopped. 
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After an interval, the timing disc actuates the device in such a way that it releases 
the loading rod and at the same time operates the relays which restart the drum and 
switch on all the lights. The timing disc and the arms which make contact with 
it are so arranged that the times at which unloading and reloading take place can © 
be adjusted at will, and the process repeated a number of times if desired. 
The wire is surrounded by a thermostat and the load applied through an 
Andrade—Chalmers constant stress bar. 


§3. MECHANICAL AND OPTICAL DETAILS 


The recording apparatus is mounted on a steel base plate 60x 45cm. The 
Muirhead 220v. phonic motor, controlled by a calibrated 50 c/s. fork, drives, 
through a 50,000 to 1 reduction gear, (1) the vertical shaft which carries the light-’ 
alloy recording drum D, 60cm. high and 21cm. in diameter, the photographic 
paper which is wrapped around it being shielded by a housing H pierced by a 
vertical slit to admit the operating beams of light, and (2) the horizontal shaft 
which carries the ebonite timing disc around which are disposed platinum studs. 
‘These studs make contact with brushes carried on two arms which can be rotated 
about an axis in line with that of the disc and can be clamped in positions read — 
off on a circular scale taken from an old spectrograph. Of the contacts with these © 
two arms, one serves to remove the load and the other to impose it again, by means 
of a device described later. 

A plan of the optical system is shown in Figure 1. The periodic flash from 
the lamp A, which records the time scale on the paper, takes place every 16 seconds, 
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a period dictated by the gears available. ‘The light illuminates the slit T, an image 
of which is thrown on the drum through the intermediacy of a cylindrical lens Ly, — 
a converging Jens L, and a front-silvered mirror M,. The horizontal length _ 
scale, providing lines spaced at about 3mm. on the drum, is produced by an 
‘mage of the scale S, reflected in the front-silvered plane mirror M,. This scale is 
made up of lines 0-5 mm. apart, so that the magnification is 6. 

To produce the trace that gives the length of the flowing wire a good point 
source attached to the lower end of the wire is required. This is provided by a 
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small steel sphere of diameter 0-4mm. attached to a strip of mirror M3, the 
| movement of which is restrained to a vertical plane by means of guide rods 
4 arranged so as to giveaminimumof friction. This mirror, being arranged to throw 
| light from the source B, rendered parallel by the lens L3, well away from the axis 


OO’, ensures a good dark background for the light proceeding to the drum from 
the minute image of the light source formed by the sphere. An aperture F 


+ confines the light to the mirror and so avoids scattered light. ‘This arrangement 


gives a fine trace which, with the magnification used, enables the length of wire 


|) to be read to within 1/50 mm. if necessary. A total extension of the wire of 10 cm. 


can be recorded. ‘To make examination of the resulting graphs easier, every 

tenth horizontal and vertical line is automatically made darker than the others. 
The systems A to L,, B to L3, C to S and L, to M, are mounted on four optical 

benches which are attached to the girder system which carries the thermostat, 


} the constant stress bar, the unloading device and associated equipment, and the 


support for the upper end of the wire. Ilford document paper has been found 
to be very suitable for the recording. Part of arecord obtained from the apparatus 
is shown in Figure 5 (Plate). 

The unloading device, shown in Figure 2, consists of two jaws J, J attached to 


Figure 2. 


| blocks B, B which are threaded to take the screws S, S. These screws are right- 
{ handed and left-handed at the opposite ends respectively, so that rotation in one 
' direction brings the blocks together and in the other direction takes them apart. 
| ‘The cross section of a jaw is shown at CC’: it is clear that a horizontal pressure on 
? the edge E will cause the jaw to rotate about the red P, the rotation being checked 


by the projection shown above the jaw. ‘Thus, when the two jaws come together, 


t they grasp thevertical rod R, which applies the load to the wire, and raise it slightly 
4 (about 0-5mm.). When they are driven apart, the rod will be gently released 
' and the load applied. 


§4. DETAILS OF ELECTRICAL OPERATION 


The complete electrical circuit diagram is given in Figure 4, the diagram 
showing the units involved and the cabling between them. 
The facilities which the apparatus provides may best be discussed by following 


‘| through the processes involved, enlarging, where necessary, on those points where 
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alternatives are offered. The light supply is 230 v. a.c., the fork and motor supply 
is 230 v. D.c., the relays and clutches 130. p.c. and certain power relays and the. 
timing lamp 24v. p.c. All supplies are controlled in the first place by the power 
relay AP (Figure 4) which must be manually closed before any power can be put 
on to the equipment, the relay holding itself from the derived 24v. supply. 
Should, for any reason, the power fail, this relay will release and will not remake. 
when the power is re-established. 

The 130v. p.c. is provided by a bridge metal rectifier and a Variac transformer, 
this supply being common to relay and timing units and to the clutches. The 
24. D.C. is provided through a mains rectifier set, the input being derived through 
the AP contacts. With AP closed and the fork set in vibration, the phonic motor 
may be started, and with the Variac output set at about 180v., both clutches will 
operate. The closing, every 16 seconds, of the cam-operated contacts CM 
operates relay BS, which releases BP: this operates the large mercury relay MR, 
which closes the flash-lamp circuit. The cam-operated circuit also contains a 
1,000-ohm counter which counts, when necessary, the number of pulses delivered. 
This circuit is used (a) to calibrate the timing against a standard clock, and (6) 
to record the duration of the ‘‘off-load”’ period, so that the accuracy with which 
the hands can be set against the scale is not a limiting factor. The exact time of 
creep under load can, of course, be read direct from the photographic record. 
In addition, the pulses delivered by the cam-operated contacts operate a P.O.-type 
uniselector which provides a higher voltage and a brighter flash at every tenth 
pulse by shorting out a resistance in the lamp circuit. 

What may be called element 1 of the relay unit is shown in Figure 3 (a). 
The relay A is operated by the contact on the arm M coming into contact with the 
stud embedded in the timing disc. This stud is earthed, and the A contact 
operates a power relay CP (see Figure 4) which drives the unloading motor forward}. 
the two jaws close, lifting the load and at the same time breaking all the light 
circuits and the circuit through the recording drum clutch. The stud is eventually 
brought by the phonic motor into contact on the second, ‘“‘reload’’, arm N, 
causing B and BB to operate the relay DP, which reverses the polarity of the field 
windings in the loading motor so that it releases the load, continuing to drive 
back until its circuit is broken by another mechanical micro-switch. At the same 
time the lighting circuits are remade. 

On the release of B, if only one interruption is to occur, the circuit is cut off 
and any further contacting, in rotation, has no effect. This is accomplished by 
a relay RA which operates when B and BB operate, and holds by one of its own 
contacts (see Figure 3(b)). When B and BB release, the B contact falls back on 
to its rest position, but as RA has not had time to release, a pulse is provided through 
RB which holds through its own contact, preventing any further operation of che 
arm circuits. It may happen that, because of the geometry of the timing disc 
arrangement, the “reload” arm N is set in advance of the “unload” arm M.° Asa 
consequence, the operation of B and BB, and the cutting out of the circuits by 
RB, would mean that unloading would never occur. Hence R is introduced 
(see Figure 3(a)) into the circuit, with the function of preventing the “load” 
circuit from operating unless A has operated and unloading has occurred.. The 


effect is that the disc will drive past the load arm N to the unload arm M, and only 
operate ‘“‘reload”’ on its next revolution. 


} i | } 
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Part of a record. showing the first few minutes of a creep curve, 
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Repetition of the Unloading Operation 


If more than one unloading operation is required, then when the load has been 
replaced after the first operation, the timing system should return to the condition 
of the outset of the experiment, with the proviso that it retains a “memory” of 
one operation. As a consequence, if S in Figure 3(c) is switched to any of the 
positions 2 to 6 (representing the number of unloading events required) the 
operation of F, which is introduced in series with RB and operates whenever RB 
does (but has no influence in the single interruption arrangement), breaks the 
circuit through the timing clutch disc. This clutch is pulled out by springs— 
it is a friction cone type of drive—and the disc spins back under the force of the 


to Arm 
Contacts 


to earth 
Pp by contact, = 
on disc 


Figure 3. 


spring mounted coaxially with the disc. A contact guided by a screw thread 
in the edge of the disc, in which it rides, moves back a conducting arm until the 
original, zero, position is reached, when the contact spindle snaps into a spring 
(“Rim Catch ’”) at the edge of the disc, the clutch circuit is restored, and a pulse 
is delivered through the contact of switch S. If, for example, the switch S is set 
at 5, K operates, in turn operating N, which holds itself and steps the feed on to 
the next relay from K, i.e. J, ready for the next pulse. Simultaneously the pulse 
operates P, whose contacts break the power to all the relays concerned with the 
primary operations of loading and unloading, so that, as far as these are concerned, 
the circuits are as they were at the beginning of the experiment. The ‘‘ memory”’ 
has been achieved by the operation of N, which remains held for the rest of the 
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experiment. This whole sequence is therefore repeated, che relays O, Q operating 
in turn until the feed is led direct to the relay H. This terminates the series, for 
the next pulse operates H direct, which, it will be seen, operates T. T holds itself 
through one of its contacts, the other contact short-circuiting the relay F, so that 
the next time an unloading operation is completed, only RB will operate, as at the 
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end of the single-interruption case; however long the disc continues to be driven, 
no further interruptions will take place. Hence a selectivity of 0-6 interruptions 
is at our disposal, with a choice of the settings of the hands anywhere within the 
360°, the ‘“‘off-load”’ period being variable, also, up to 360°, which corresponds 
to one revolution of the recording drum. 
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Additions of ‘‘ Whole-revolutions’’ to the Off-load Period 


The off-load period has been made automatically controllable for a period 
longer than that represented by one revolution (i.e., approximately 14 or 3 hours, 
depending on the gear ratio used). 

This has been effected by incorporating a second contact on the timing disc, 
diametrically opposite to the other contact, i.e., at the bottom of the disc in the 
zero position shown in Figure 4. Being atagreater radial distance from the centre, 
it does not contact the setting arms, but instead, a particular contact fixed to the 
base plate. The relay system which provided this facility may be represented as 
shown (Figure 3 (d)), in which the “barring” relay R to the B and BB relays (i.e., 
R must be operated, by A, before the B and BB relays can operate) is now operated, 
not by A, but by another relay Z. 

Consider the circuit shown in Figure 3 (d) when the contact point P is connected 
to earth and then disconnected (that is, an earthed contact is driven across the arm) 
a number of times. 

Make 1. W and U operate. 

Break 1. Ureleases: Wisslowto release, sothat Y receivesa pulse and operates 

and holds through S so that both S and Y are now held. 
Make 2. Because Y is operated, this pulse operates V and W. 
Break 2. VY releases. W is slow to release, so that X receives a pulse, operates 
and holds through its own contact, switching the input from 
P on to Z. 

Make 3. Z operates, which operates R and hence makes it possible for B and 
BB to operate on their next contact of the ‘‘reload”’ arm. 

By introducing a switch into the circuits, the point at which this process starts, 
and, therefore, the number of pulses required to operate Z, may be selected. 
If a period of off-load longer than 12 hours, which is the maximum under auto- 
matic control, is required, it can be given by manual switching. 

The two setting hands and the two switches, of which one selects the number 
of off-loadings and the other adds time increments to the off-load period, provide 
the required pre-set facilities. 

When the experiment has been terminated by one complete revolution of the 
drum, a micro-switch in the drum-housing operates AS and hence releases AP, 
cutting off all the power supplies. 
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ABSTRACT. Amethod of dete1mining the stress-strain relation of materials when stresses 
are applied for times of the order of 20 microseconds is described. The apparatus employed 
was a modification of the Hopkinson pressure bar, and detonators were used to produce 
large transient stresses. Thin specimens of rubbers, plastics and metals were investigated 
and the compressions produced were as high as 20°% with the softer materials. It was found 
that whilst Perspex recovered almost as soon as the stress was removed, rubbers and poly- 
thene showed delayed recovery, and copper and lead showed irrecoverable flow. The 
phenomenon of delayed recovery is discussed in terms of the theory of mechanical relaxation 
and memory effects in the material. 


§1. INTRODUCTION 

HE mechanical behaviour of materials at high rates of loading is of interest 

from two points of view. Firstly, the engineer, who wishes to use the 

materials under conditions where they may have to withstand sudden impacts, 
requires to know how their mechanical properties depend on the rate at which 
stresses are applied. Secondly, since the variation of the stress-strain diagram 
with loading rate is related to relaxation processes taking place on a microscopic 
scale in the material, the results are of interest to the physicist who is studying 
the relation between physical problems and molecular structure. 

The present investigation was undertaken largely from this second point of 
view and the first experiments were carried out with specimens of several rubbers 
and plastics. Earlier work (Naunton and Waring 1938, Alexandrov and 
Lazurkin 1940) has shown that the mechanical behaviour of these materials is 
highly dependent on the rate of application of stress, aad it was considered of 
interest to investigate their mechanical properties when they are being submitted 
to large stresses which are maintained for only a few microseconds. Some experi- 
ments on the behaviour of copper and lead under these conditions have also been 
carried out. 

Two problems are encountered in such an experimental investigation of the 
mechanical behaviour of maierials at very high rates of loading. These are 
associated with inertia effects in the apparatus and with the recording of transient 
stresses and strains. 

As the rate of loading is increased, the acceleration of any moving parts of the 
straining apparatus begins to require forces comparable with those necessary to 
deform the specimen. It is often impossible to separate these inertia effects 
from the effects due to the physical properties of the material under investigation, 
and at the highest rates of loading the inertia of the specimen itself will result in 
a non-uniform distribution of stress along its length. In the work described here 
these difficulties were overcome by the use of very thin specimens and the stresses 
were applied by a simple mechanical system for which the inertia effects could be 
calculated in terms of the propagation of stress waves. 
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The natural period of ordinary stress gauges is far too long for them to respond 


‘to forces which are applied for only a few microseconds, and special techniques 


must be used for recording such transient stresses. Hopkinson (1914) was the 
ficst to devise a satisfactory method of doing this and the technique was later 


developed by Landon and Quinney (1923). More recently R. M. Davies (1948) | 


has made a critical study of the method and has described a modified form of the 


apparatus in which the measurements are made electrically. This overcomes 


many of the disadvantages associated with the earlier experimental work. 

Hopkinson’s original apparatus, which has become known as the Hopkinson 
pressure bar, consisted of a cylindrical steel bar suspended ballistically. The 
pressure to be measured was applied to one end of the bar whilst at the other end 
a cylindrical pellet, known as the time piece, was wrung on. ‘These pellezs were 
of the same diameter as the bar and of the same type of steel. The end of the bar 
and one surface of the pellet were ground flat so that the pellet could be attached 
to the end of the bar with very little grease. 

The principle of the method is that when a pulse of compression travels down 
the bar it is transmitted through the greased joint without change of form. When 
it reaches the other face of the pellet, which is a free boundary, it is reflected as a 
pulse of tension. Since the greased joint is unable to withstand any appreciable 
tensile stress, the pellet flies off as soon as the reflected pulse builds up any tension 
across the interface. The momentum trapped in the pellet corresponds to a 
section of pulse the length of which is twice that of the pellet. If the experiment 
is repeated with pellets of different length, and the momentum measured in each 


case by capturing the pellet in a ballistic pendulum, the nature of the pressure- 
‘time relation for the pulse can be investigated. 


Although Hopkinson’s original method has the advantage of simplicity and 
has been used with success by the Research Department, Woolwich (Robertson 
1921), for measuring the pressures set up by the detonation of various explosives, 
it suffers from several serious limitations. Firstly, it gives only a series of pressure— 
time integrals for different sections of the pulse. From these the maximum 
amplitude of the pressure, and the time for which any value of the pressure was 


-exceeded may be deduced, but the actual pressure—-time curve is not determined. 


Secondly, the force necessary to break the greased joint introduces an unknown 
variable into the experiments, and this precludes the use of the apparatus for the 
measurement of pulses of small amplitude. Finally, the method assumes that 
the pulse is propagated down the bar without appreciable change of form: this 
is true only for pulses which are long in comparison with the diameter of the bar 


-and within which there are no sudden changes in pressure. 


In the Davies method (Davies 1948) the first two of these difficulties are 
overcome by eliminating the pellet and greased joint and obtaining a continuous 
record of the very small displacement of the free end of the bar. This is done by 
making the free end of the bar the earthed conductor of a parallel-plate condenser 
which operates as a condenser microphone. The amplified output from this 
microphone is fed into a cathode-ray oscillograph and a photographic record made. 
From this record the displacement-time curve may be obtained and if the applied 
pressures do not exceed the elastic limit of the steel the pressure-time curve can 
be deduced. Davies has also investigated the theory of the propagation of pulses 
down cylindrical bars and shown under what conditions the distortions become 
too severe for the method to be used satisfactorily. 
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In the work described here an apparatus based on that designed by Davies has. 
been used for making stress-strain measurements. The specimens, which were. 
in the form of thin discs, were placed between the flat faces of two cylindrical 
steel bars. The transient pressure was applied by firing a detonator at the end 
of one of the bars and the displacement at the free end of the other bar was measured 
with a parallel-plate condenser microphone. From this the pressure on the. 
specimen could be calculated. A cylindrical condenser microphone was also 
fitted around the bar between the detonator and the specimen; this was used to 
measure the amplitude of the pressure pulse arriving at the specimen, and the- 
deformation of the specimen could then be deduced. 


§2. APPARATUS 
(1) General Description 
Figure 1 shows the general arrangement of the apparatus; this is similar to 


that described by Davies except that the bar is in two parts and a second condenser” 
microphone and amplifier are introduced. 


; Cylindrical Parallel Plate 
é Inertia Condenser Condenser 
Aawil Collar Switch Microphone Microphone 
/ M 


Detonator Collar Specimen 


Cathode 
Ray 
Oscillograph 


Figure 1. General arrangement of apparatus. 


(ii) The Pressure Bar 


The bar and extensions were of silver steel (yield point 40 tons/in?) 1 in. in. 
diameter. ‘The bar was 6 feet in length and extensions 4 in., 6 in. or 8 in. in length 
could be attached with a brass collar, the specimens being lubricated and placed. 
between the flat faces of the bar inside the collar. 'The ends of both the main 
bar and the extensions were carefully ground and the collar was made to be a 
good sliding fit over both bars. In order to protect the pressure end of the bar, 
the detonators were fired against hardened steel anvils of the same diameter as. 
the bar. ‘These anvils, which were replaced frequently, were cylinders 14 in. 
in length with their ends ground flat. A close-fitting steel collar held them against: 
the end cf the bar, a very small quantity of grease or heavy oil being used between 
the flat surfaces to ensure uniform contact. 

The detonators, which could be fired electrically, were specially prepared by 
the Research Department of the Nobel Division of Imperial Chemical Industries 
Limited, and consisted of flat-based aluminium tubes 0-25 in. in diameter containing: 
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0:5 gm. of tetryl loaded and pressed in increments of 0-1 gm., and primed with 
0:35 gm. of a composition of lead azide and lead styphnate. The results obtained 
with these detonators were found to be consistent to one or two per cent. 

The detonators were fitted to the anvils by means of cylindrical polythene 
holders down the axes of which holes of suitable diameter had been drilled. The 
holders were attached to the anvils with adhesive tape so that the flat base of the 
detonator was firmly in contact with the face of the anvil. The bar was suspended 
ballistically and the momentum communicated on firing the detonator was 
measured with a light metal pointer which gave the maximum amplitude of swing. 
The various electrical connections were made with thin flexible wire so that they 
did not hinder the motion of the bar. 


(ii) The Condenser Microphones 


Figures 2 and 3 show the construction of the two condenser microphones 
used with the apparatus. The cylindrical microphone shown in Figure 2 consists 


Figure 2. Cylindrical condenser microphone. 


Figure 3. Parallel-plate condenser microphone. 


of a circular disc of mild steel A din. thick, through which a concentric hole 
1-020 in. in diameter has been drilled; this disc is held between two mild steel 
collars B by four bolts. ‘Two Perspex washers C were used to insulate the disc 
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from the collars as this material was found to give good electrical insulation and 
mechanical rigidity. ‘The collars were machined to slide easily over the bar D, 
the bore being slotted at four places P to allow feeler gauges to be inserted. The 
microphone was placed on the bar, the insulated ring was set concentrically by 
means of feeler gauges and the bolts were then tightened. This left an annular 
gap of 0:010in. between the bar and the ring. 

This condenser microphone was normally placed 7in. from the end of the 
main bar and the insulated disc was charged to two or three hundred volts through 
a high resistance. The condenser feed-unit is shown in Figure 4. Whilst the 
compression pulse from the detonator was passing through the microphone, a 
lateral expansion of the bar took place, which caused a change in the capacity of 
the cylindrical condenser consisting of the bar and ring. ‘his resulted in a 
small change in voltage across the condenser which was amplified and fed to one 
beam of a double-beam oscillograph as shown in Figure 1. 

The parallel-plate condenser microphone shown in Figure 3 is similar to the 
one described by Davies (1948), except that provision has been made to set the 
insulated plate accurately by adjusting the bolts. The insulated metal plate 


R, c, 
l l 
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Figure 4. Condenser feed unit. 
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A is held in the cylindrical block of polythene B which is bolted to a mild steel 
collar C, a soft rubber disc D being placed between the polythene and the collar. 
The collar was made to slide easily over the ends of either the extension bar E 
or the main bar. ‘The separation between the plate and the bar was adjusted 
by inserting ring gauges between the collar C and the small ridge F. The micro- 
phone was connected through a feed unit similar to that shown in Figure 4. The 
output was then amplified and fed to the other beam of the oscillograph (see 
Figure 1). 

For both types of microphone the movements in the bar are converted into 
voltages which are proportional to the corresponding displacements so long as 
these displacements are small compared with the distance between the plates of 
the condenser units. For larger displacements the relation ceases to be linear. 

As can be seen in Figure 4, the condenser microphone C, is connected across 
a much larger capacity C, and these are both charged through a high resistance 
R,, Cz being a large isolating condenser between the feed unit and the amplifier. 
The value of R, is sufficiently large for the charge across the two condensers in 
parallel to remain constant whilst the value of C, is being changed by the pressure 
pulse. ‘T’hus, if we write C for C\+C,, and c, is the change in capacity of Cj, 
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V the original potential difference across the condensers and w the change in 
potential difference, we have 

O=VC=(V—a(Ereg) gh engines (1) 
and if c, is small compared with C (as it was with both condenser units), we may 
write 


OCG Ee MPP iia te AL es (2) 
For the parallel-plate condenser unit the capacity C, is given by 
OF NORE) ct ONS 6 Wire nee Po: amare se (3) 


where A is the area of the plates, D the initial distance between them and ¢ a 
correction term to allow for edge effects. This term is small when D is small 
compared with the radius of the plates. If wetake C, = B/D where B= A(1+4)/47,. 
B will remain appreciably constant for small changes in D. 

If the end face of the bar moves a distance x to produce the change in capacity 
C,, we have C, +c, =B/(D—x), so that 


(5 dad oh COLES) D ee ie emis ess (4) 
or substituting in (2) 
: = VBx(D—x)CD=VCw(D—2)C. = 11s (5) 
If x is very small compared with the separation D 
VG SADR ON Tare oy Med eee (6) 


so that the voltage is proportional to the displacement. 

Since the annular gap in the cylindrical microphone is small compared with 
the radius of the bar 7, the capacity of the condenser is here again inversely 
‘proportional to the width of this gap and equation (5) may be used. _D is in this 
case the width of the gap, x the change in radius of the bar and, if edge effects 
are neglected, B=rl/2 where / is the thickness of the insulated disc. 


(iv) Amplifiers 

The output obtained from both types of condenser microphone described 
above will be proportional to the voltage applied across them, but even if this is 
increased to 1,000 v. the output resulting from the small movements of the bar 
produced by the pressure pulse will still be only a small fraction of a volt. In 
order to use this to produce a measurable deflection of the beam of a cathode-ray 
oscillograph a high-gain amplifier is therefore necessary. ‘The design of such 
an amplifier which will give undistorted amplification of short transient pulses 
is a matter of some difficulty since stray capacities result in poor high-frequency 
response and attempts to compensate for this by using partially inductive loads 
in the anode circuits may easily lead to the setting up of parasitic oscillations. 
Fortunately, satisfactory wide-band amplifiers of high gain have been designed 
recently (Valley and Wallman 1948) and one of these, the No. 234 Pulse Amplifier, 
developed by the Telecommunications Research Establishment, was used in the 
present investigation. ‘This amplifier gives an even voltage amplification from 
audio frequencies up to frequencies of several megacycles per second and has a 
maximum gain of about 30,000. 

The amplifier is designed to work only with pulses of one sign but this did not 
matter with the parallel-plate microphone as the displacement of the bar resulting 
from a pressure pulse is always in one direction. With the cylindrical unit, 
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oscillations in the displacement are found to occur and in order to observe these 
a wide-band amplifier of smaller gain, which did not suffer from this limitation, 
was employed. 

The connections between the feed units and the amplifiers were made with 
low-capacity coaxial cable and this was also used for connecting the outputs of 
che amplifiers to the Y plates of the double-beam oscillograph. 


(v) Cathode-ray Oscillograph 
A Cossor oscillograph Unit Model 3402, employing a 3259J double-beam 


tube, was used in these experiments. This tube can be ru at a maximum anode 
voltage of 3,000 v. and was found to give an adequate writing speed for the experi- 
ments when fast photographic plates and a large aperture camera lens were 
employed. 

The output from each amplifier was connected to the Y deflector plates of 
one of the beams, the other deflector plate in each case being at earth potential. 
Provision was made to substitute an oscillator of accurately known frequency 
on to one of the Y plates, and this was used as a time calibration for the records. 

The X plates were connected across a balanced time-sweep unit which, when 
triggered, resulted in a single traverse of the two beams across the screen. At the 
same time it applied a transient positive potential to the modulating grid of the 
ccathode-ray tube, which brought the beam to maximum brightness during the 
sweep. 

The camera lens used had an aperture of f/2-9 and was of 8 in. focal length, 
a magnification of 0-5 was generally employed. In order to ensure that no 
distortion was being introduced by the camera, a graduated grid was placed in 
the plane of the cathode-ray tube screen and photographed; the plate was then 
measured with a travelling microscope. Ilford H.P.3 quarter plates were used 
and these were found to give a considerably higher writing speed than other plates 
with the J screen tube. The plates were measured with a Beck travelling micro- 
-scope which gave traverses in two directions at right angles. 

In carrying out the experiments the brilliancy control of the cathode-ray 
oscillograph was adjusted so that the spots just disappeared, the room was then 
placed in darkness and the camera shutter opened. When the detonator was 
fired it actuated an inertia switch on the bar which triggered the sweep unit giving 
a single traverse of the two beams. 


(vi) Sweep Unit 

The circuit of the sweep unit, together with the values of the resistances and 
condensers normally used, is shown in Figure 5. The sweep unit is designed to 
give a large transient voltage between X1 and X2 when a positive pulse is applied 
toterminal T. ‘To prevent defocusing of the beam this voltage is balanced about 
earth potential. As soon as the unit is triggered it also applies a positive voltage 
to the modulating grid of the cathode-ray tube from terminal M, and this voltage 
is maintained positive for the duration of the sweep; it then falls back through 
zero to a small negative value and then asymptotically back to zero again. V1 is 
a rectifying valve and V2 a gas-filled or mercury vapour relay valve. For sweep 
voltages up to 500 v. a GT1C Osram valve with indirectly heated cathode was 
used, ‘I'o obtain larger sweeps a BT19 mercury vapour thyratron working at 
1,000 v. was employed. 
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The potentiometer R6 could be set so that the sweep unit triggered itself 
and gave a repeating stroke. It was used in this position whilst the other controls 
were being adjusted. 

In order to apply a triggering voltage to terminal T some microseconds before 
the pressure pulse reached the cylindrical microphone, an inertia switch was 
fitted to the bar 6in. from the microphone (see Figure 1). Several types of inertia 
switch were tried in attempts to effect triggering at a reproducible time before the 
head of the pressure pulse reached the specimen. With all mechanical types an 
uncertainty of 2 or 3 microseconds was found to occur and a simple ring switch 
similar to that described by Davies was finally employed. This is shown dia- 
grammatically in Figure 1 and consists of a metal tube pressed round a plastic 
insulating collar which slides easily on the bar. 


C9. RI4 
eT XI 
R10 
a Seal C7 Ri2 
230V 
oe; 
4 RIB 
X2 
cl0 RIS 
Figure 5. Sweep unit. 
R1, 14,15 30,0009. C1) 0-5 EF. 
R2, 3 0:5 MQ. C2 10 HF. 
R4, 5 2MqQ. C3 50 HE. 
R6, 7, 8, 9 100,000 Q. C4, 6, 7, 
R10, 11 20,000 Q. S910 0-01 pF. 


RI2ZA13 1 MQ. C5 0-001 pr. 


It was found that the flash from the explosion of the detonator could be used 
to trigger the sweep unit satisfactorily by letting the light fall on a photoelectric 
cell with a one-valve amplifier. By this means results reproducible to within 
one microsecond were obtainable and this method was employed for some of the 
experiments. 

Figure 6 shows a photograph of the traces obtained in an experiment when the 
bar and extension were wrung together with a thin layer of grease between them. 
The top and middle traces are the amplified outputs from the cylindrical and 
parallel-plate microphones respectively. ‘The bottom trace is a timing wave of 


5-30 psec. period. 
(vii) Static Stress—Strain Apparatus 


Figure 7 shows a diagrammatic representation of the apparatus used to obtain 
stress-strain curves for the materials under static loading conditions. ‘These 
were carried out for comparison with the results obtained for the Hopkinson bar 
apparatus. ‘The specimens were placed in the hydraulic ram between the flat 
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end of the piston which was 1 in. in diameter and a stop of the same diameter, the: 
faces of the specimen being lubricated to allow free lateral movement. ‘The stress 
was measured with the pressure gauge which gave the oil pressure on the piston. 
The change in thickness of the specimen was measured with a travelling micro- 
scope which read directly to0-01 mm. The pressure was built up with the pump, 
the booster being used as a fine control. 


§3. THEORY OF THE METHOD 


If a plane compression pulse is propagated along a cylindrical bar and the 
pressure does not exceed the elastic limit of the material, the velocity of propagation 
will be c where c=1/(E/p), E being Young’s modulus for the material and p its 
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Figure 6. Cathode-ray oscillograph traces showing outputs from the two condenser microphones: 
and 5-3usec timing wave. 


Figure 7, Static stress-strain apparatus. 


density. Also if the particle velocity at any point is w and the amplitude of the 
pressure 1s p, we have p=pcu. (‘These relations assume that the pulse is pro- 
pagated without distortion, the pressure being uniform over any given cross 
section; this is true only for pulses which are long in duration compared with the 
time taken for the pulse to travel a distance equal to the radius of the bar.) 

At the free end of the bar the compression pulse will be reflected asa pulse of 
tension. ‘T’his is of the same form as the outgoing pressure pulse, and the velocity 
of the end of the bar is therefore twice the corresponding particle velocity associated 
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with the pulse travelling along the bar. Thus, if we denote the displacement of 
the end of the bar, at time ¢, by € we have 


d. 
p=decS. eet (7) 


The output from the parallel-plate microphone gives the relation between 
€ and ¢ and by differentiating this curve the pressure-time relation for the pulse 
may be derived. 

The radial displacement of the bar ¢ is given by 


Gas OT PL a en eee aay (8) 


where o is Poisson’s ratio and r is the radius of the bar. Thus the cylindrical 
condenser microphone gives the pressure-time curve directly, when the value of 
or/E for the bar is known. 

Davies (1948) has considered the two ways in which distortion appears when 
the pulse is of a length comparable with the cross section of the bar. He has done 
this in terms of the elastic equations given by Pochhammer (1876), Chree (1889) 
and Love (1927) for the propagation of longitudinal vibrations of a circular cylinder. 
This treatment shows that for such short pulses considerable dispersion takes place, 
the Fourier components of short wavelength travelling with lower velocities than 
those of longer wavelength. ‘This results in a lengthening of the pulse as it 
travels down the bar. The other form that the distortion takes is in a non- 
uniformity of the distribution of the pressure and displacement across any section 
of the bar. When measurements are made with the parallel-plate microphone 
this results in an amplification of the high-frequency components. 

Davies estimates that when a bar 2ft: long and 1 in. in diameter is used for 
measuring pulses of 20 sec. duration the distortion causes errors of between 
2 and 3% in the values obtained. He also describes experiments which confirm 
that this is the order of the error observed. he distortion results in larger errors 
when a cylindrical condenser microphone is used and Davies suggests that this 
method can only be used profitably with pulses which are very long compared 
with the diameter of the bar. 

In the method described here for obtaininy stress-strain measurements, any 
distortion introduced by the main bar does not affect the validity of the results, 
but only changes the form of the stress cycle to which the specimen is subjected. 

If a sufficiently thin specimen is inserted in the split Hopkinson bar, described 
in the previous section, the pressure 1s effectively the same throughout the specimen 
during the passage of the compression pulse, and this will be true so long as the 
thickness of the specimen is small compared with the wavelength of the shortest 
operative wave in the Fourier spectrum of the pulse. This pressure is com- 
municated along the extension bar to the parallel-plate condenser microphone. 
On differentiating the displacement-time curve, which is obtained from this 
microphone, the stress-time relation for the specimen may be derived using 
equation (7). 

If the displacement-time relation for the initial pulse travelling down the 
main bar is also known, the strain-time relation may be determined and the 
stress-strain curve can then be derived. Let the pressure of the incident pulse 
at any time f be ,, of the transmitted pressure pulse which acts on the specimen. 
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py, and of the reflected pulse, which is a pulse of tension travelling back along the 
bar —pp, then 


Pie Pere CP ee (9) 


If we denote the displacement of the face of the main bar in contact with the 
specimen at time ¢ by €, we have 


ie 
é,= sq |, (b+ nda 1 oa (10) 


and if the displacement of the face of the extension bar in contact with the 
specimen is £5, 


1 ct 
a= =| pra wae (11) 


Now if the thickness of the specimen is s, the fractional strain s is given by 
(€, — &)/z, hence from (9), (10) and (11) we have 


yy t 
s= = | (or pat Mie (12) 
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is the displacement measured by the parallel-plate microphone with the specimen 
between the bar and the extension, and 


2 ft 
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the displacement produced by the pulse when no specimen is present, the two 
bars being in contact. Thus 
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so that if we know the displacement-time relations for & and €» we can obtain 
the strain—time curve and, since the gradient of the &, curve gives us the stress— 
time curve, we can derive the stress-strain relation for the material under the 
conditions of the experiment. 

As mentioned above, the distortion introduced by the bar causes the 
values of the pressure amplitude obtained with the cylindrical microphone to be 
rather unreliable. Nevertheless, the experimental results showed that the 
pressure-time integral, which is all that is required to determine the strain in 
the specimen, agreed closely with the values obtained directly with the parallel- 
plate microphone. 

In the stress-strain determinations described here, however, it was generally 
found most convenient to determine é, and €) in separate experiments with similar 
detonators, the two bars being wrung together with a little grease in the first case, 
and the specimen lubricated and inserted between the bars in the second case. 
The cylindrical-condenser microphone, did, however, serve two useful purposes. 
Firstly, it gave a convenient method of fixing the same zero in time for the two 
experiments, since measurements along the time scale could be made from the first 
intercept of the cylindrical microphone trace on the horizontal axis (Figure 6, 
top trace). ‘This could be done with considerable accuracy whereas, since the 
displacement begins rising very slowly (Figure 6, middle trace), the arrival of the 
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pressure pulse cannot be fixed at all accurately. Secondly, the trace from the 
cylindrical microphone gave an accurate check that the pressure pulses obtained 
with different detonators were of the same shape and amplitude. 

The forward acceleration of the specimen is proportional to the difference in 
pressure between the two faces and, so long as the specimen is sufficiently thin, 
the inertia of the specimen may be neglected from this point of view in determining 
its stress-strain behaviour. The radial motion, however, depends on the time 
rate of change of the longitudinal strain and a correction must be applied for the 
radial kinetic energy however thin the specimen may be. 

If we take p, as the stress required to produce a fractional strain s when no 
kinetic energy is given to the specimen, p,, to be the measured pressure with the 
specimen in the split Hopkinson bar apparatus, and G to be the radial kinetic 
energy, we have 


Az | Ca ae ce aie ey OM ee (14) 


where A is the cross-sectional area and z is the thickness of the specimen. 
Differentiating with respect to s we obtain 


dG ; 
Peas rae oieueueiele (15) 


Now if the specimen is expanding uniformly and friction between it and the bars 
is neglected, the radial kinetic energy may be obtained by integration and is 
given by 

(Chea Von 1 Ga alanle aed al LU’ ge fe (16) 
where p, is the density of the specimen and I, is the radial velocity of a point on 


the circumference. 
If o, is Poisson’s ratio for the specimen and , its radius, we have 


V.zorS, 

and substituting in (16) and differentiating we obtain 

dG d*s 
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hence from (15) 

pega” ‘ 
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This correction was found to be not more than a few per cent in most of the 
experiments and is only important when the rate of straining is changing rapidly, 
which, in these experiments, was at the beginning of the stress cycle. Nevertheless, 
having to use this correction does introduce a certain amount of error into the 
results, in that the value of o, at these very high rates of loading is not known 
accurately, and the assumption that the expansion is uniform with no frictional 
effects is not necessarily justified. ‘The double differentiation of the strain-time 
curve may also result in appreciable errors when the correction term is comparable 
with the actual applied stress. 

The value of the correction term is proportional to the square of the radius of 
the specimen and if the method is to be used at higher rates of straining, smaller 
specimens would have to be employed. 
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§4. CALIBRATION OF THE APPARATUS 


The foregoing treatment to the experimental method assumes that the 
apparatus fulfils the following conditions : (a) the extension bar does not introduce 
any further distortion to the stress pulses transmitted through it; (6) the pulse is 
transmitted across a greased joint between the bars without distortion; (c) the 
brass collar connecting the two bars does not transmit any appreciable stress- 
wave from the main bar to the extension bar. 

A series of experiments were carried out to test the extent to which these 
conditions were fulfilled by the apparatus, ‘To investigate conditions (a) and (d) 
the collar and extension bar were removed and the parallel-plate microphone was _ 
set up at the end of the main bar. Records were then taken of the pulses from 
detonators, and the results compared with those obtained when the microphone 
was at the end of an extension bar wrung on to the main bar with a little grease. 
It was found that the results obtained with 4in., 6in. and 8in. extension bars all 
gave the same displacement-time curve as that found when the microphone was 
at the end of the main bar. 

Thus conditions (a) and (b) were fulfilled for the pressure pulse which had 
arrived at the end of the main bar. The method also assumes, however, that the 
modified pulse, produced when a specimen of some other material is inserted — 
between the bars, travels down the extension bar without distortion. In the 
experiments it was found that these pulses were longer than the original pulse 
from the detonator and that the changes in pressure were generally more gradual. 
From Davies’ treatment one would therefore expect that these pulses would travel 
down the extension bars with even less distortion than the pulse from the detonator.. 
This could not be tested directly, but it was found that here again identical 
results were obtained whichever length of extension bar was employed. 

It is perhaps worth mentioning that the results obtained when the microphone 
was used with extension bars differed in two minor respects from those obtained 
with the microphone at the end of the main bar. The middle trace of Figure 6. 
shows the trace from the parallel-plate microphone when it was placed at the end of 
a 6in. extension wrung on to the main bar. It may be seen that the displacement 
rises to an approximately constant value at which it stays for about 35 psec., 
after which it begins to rise again. 

This second rise results from’ the reflected pulse travelling back along the 
extension bar, being reflected at the greased joint as a pulse of compression, and 
finally arriving at the condenser microphone as a pulse similar in form to the 
initial one. (The whole of the second pulse was not obtained in the experiment 
shown in Figure 6, as the amplifier was overloaded before the maximum dis- 
placement had been reached.) 

When the microphone was placed at the end of the main bar no such second 
rise could be observed since it would arrive after the pulse had twice travelled the 
full 6 ft. leugth of the bar; it would thus be well outside the duration of the trace. 
The other difference between the traces obtained when the microphone was at 
the end of the main bar, and when an extension bar was used, is that in the former 
case instead of the trace being approximately flat after the main pulse had passed, 
it showed a series of oscillations. Such oscillations have been observed by Davies 
and shown by him to be due to the distortion introduced by the bar; this dis- 
tortion results in a lagging behind of the high-frequency components of the pulse. 
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‘These oscillations may also be seen in the top trace of Figure 6 which corresponds 
to the radial displacements of the surface of the main bar when the pulse is passing. 
‘The theoretical treatment shows that these will be more marked than the cor- 
responding longitudinal vibrations in the bar and it can be seen from the Figure 
that in this case the vibrations are comparable with the main pulse. 

The reason these oscillations do not appear when an extension bar is used is 


that the greased joint does not transmit tensile stresses. The bars consequently 


separate slightly as soon as the first oscillation in the stress pulse builds up an 
appreciable tension across the greased joint; this joint thus acts as a mechanical 
rectifier. 

To test whether condition (c) was fulfilled, the two bars were separated in the 
collar by an air gap of 0-020 inches (this gap being several times larger than the 
calculated maximum displacement of the end of the bar produced by the pulse). 
‘The amplifier of the parallel-plate microphone was then set at maximum gain and 


a polarizing voltage of 360 v. was used across the microphone. ‘The sensitivity 


of the microphone was then about 30 times as great as it was in the stress-strain 
determinations. A detonator was fired at the end of the bar and the cathode-ray 
oscillograph traces photographed. ‘The trace from the parallel-plate condenser 
microphone was perfectly flat, showing that no appreciable motion had been 
communicated to the extension bar through the brass collar. 

Before the apparatus could be used for stress-strain measurements the relation 
between the values of the deflections on the cathode-ray oscillograph records and 


‘the corresponding displacements of the end of the bar had to be determined. 


In addition, as may be seen from equation (7), in order to calculate the pressure 


from the gradient of the displacement-—time curve, the value of pc for the material 


of the bar had to be ascertained. 
The relation between the voltage output of the condenser microphone system 


_and the displacement at the end of the bar is given in equation (5) and this equation 
‘may be rewritten 
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If unit deflection on the photographic plate corresponds to an output of N volts 
‘from the condenser microphone and feed unit, we have 
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where y is the deflection on the photographic plate corresponding toa displacement 
x at the end of the pressure bar. Now N may be determined if the gain of the 


amplifier, the sensitivity of the cathode-ray oscillograph and the photographic 


reduction are known. If the capacity of the parallel-plate condenser unit, C,, 


can also be determined for a known value of D, either by calculation or with a 
-capacity bridge, equation (20) may be used to calculate the displacements from the 


measurements of the cathode-ray oscillograph records. 
As this method involves the accurate determination of all the above constants, 


it was considered both simpler and more accurate to determine the quantity 


C,\V/CND=K, experimentally. From this value of K, and D as measured with 
gauges, equation (20) could be used to convert the values of y into displacements. 
If the value of y for one known displacement x can be measured the value of 
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K may be derived. Now, by integrating equation (7) it can be seen that when a 
pressure pulse is reflected at the end of the bar the maximum displacement » is. 
given by 


But if M is the momentum associated with the pulse, and A is the cross-sectional 
area of the bar, we have 


nus from (21) and (22) 


Hence, if the momentum WM associated with a pressure pulse is known, the 
maximum displacement of the end of the bar when the pulse is reflected may be 
calculated. 

Relation (23) can thus be used to determine the value of © for one known 
value of y, and K can hence be determined. ‘To do this a detonator was fired 
at the end of the bar, the maximum amplitude of swing was determined, the 
cathode-ray oscillograph traces were photographed and the gain of the amplifiers 
measured. From the amplitude of swing, the mass of the bar and its period of 
oscillation the value of M was calculated. Hence the value of x was derived using 
the values of p and c, which, as described later in this section, had to be determined 
for the pressure measurements. 

Substituting x and the corresponding value of y (measured on the cathode-ray 
oscillograph record) in equation (20), the value of K was determined. Since K 
is proportional to the gain of the amplifier, which may vary with time, in the later 
work the gain was measured in each experiment and the value of K changed 
accordingly. : 

To ensure that the momentum was entirely that of the short pressure pulse 
and that the amplitude of swing of the bar had not been increased by the general 
rush of gas following the explosion, the following experiments were carried out. 
The condenser microphone and the extension bar were removed and time pieces 
of different lengths, similar to the ones used by Hopkinson (1914) were wrung 
on to the bar. Detonators were then fired at the other end of the bar and the 
amplitude of swing measured. It was found that with a time piece 2}in. long 
there was no measurable deflection of the bar, the entire momentum having been 
trapped in the time piece. 

It may be seen from equation (20) that in order to convert the values of y into 
displacements, the value of D is required as well as the value of K. D was 
determined with mechanical gauges and there may have been an inaccuracy of 
one or two per cent in the value taken. y/D was, however, always small compared 
with K and the error introduced into the values of x by this factor was consequently 
negligible. 

As was mentioned earlier, the value of pc of the material of the bar is required 
to convert the displacement measurements into pressures. The density p of the 
steel was measured by weighing a section of known dimensions and was found to be 
7-89 gm/cm®. The velocity of propagation c of the pressure pulses was measured 
by comparing the interval of time between the arrival of the pulse at the end of 
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the main bar, and at the ends of extension bars of different lengths. Records were 
taken with the parallel-plate condenser microphone first at the end of the main 
bar and then at the end of an extension bar, the cylindrical microphone being 
kept at the same position for all the experiments. The velocity of propagation 
c could then be calculated from the difference in separation between the pulses 
from the two microphones with and without the extension bar (cf. Figure 6). 
The mean value obtained for the velocity was 5,320 metres/sec. This value was 
checked by setting a length of magnetized bar in longitudinal resonance by means 
_ of an oscillator and energizing coil, and measuring the frequency, and it was found 
that the velocity of sound measured in this way agreed closely with the velocity 
of the pulses. 


§5. EXPERIMENTAL RESULTS 
(i) General Procedure 


The specimens used for the experiments with the pressure bar were in the form 
of circular discs of uniform thickness. In order to allow for the lateral expansion 
of the specimens when they were being compressed, their initial diameters were 
made slightly less than that of the bar, the maximum permissible diameter in 
each case being determined by preliminary experiments. To reduce friction 
between the specimens and the bars both surfaces of each specimen were lubricated 
before it was placed in the apparatus. A thin layer of yrease was generally used 
except with the rubber specimens, when glycerol was used. A number of different 
oils and greases were first tried in order to ensure that the viscosity of the lubricant 
did not materially influence the results. Similar results were obtained with all 
the lubricants used; when, however, no lubricant was used, it was found that 
larger stresses were required to obtain a given strain, which indicated that there 
was then a considerable amount of friction between the specimen and the bars. 

For each group of experiments the gains of the amplifiers were measured, and 
in addition one or two experiments were included in which the main bar and the 
extension bar were wrung together with no specimen between them. The 
traces obtained under these conditions were useful as a calibration of the micro- 
phones and as a check that the apparatus as a whole was operating satisfactorily. 

The results of the ballistic calibration of the pressure bar, described in the 
previous section, showed that when the pressure pulse from one of the detonators 
used was reflected at the free-end of the main bar, this end surface moved forward 
a distance of 0:12mm. _ In the limiting case of a very soft specimen the movement 
of the end face would be very nearly as much as for a free end. The other face of 
the specimen in contact with the extension bar would remain almost fixed, and the 
maximum fractional strain would then be given by the ratio of 0-12 mm. to the 
initial thickness of the specimen (cf. equation (13)). 

With harder materials both surfaces of the specimen move and the strains are 
consequently smaller; nevertheless, by using specimens of different thicknesses, 
the rate of straining and the maximum fractional strain can be varied within 
certain limits. These limits are set on tne one hand by the fact that, for very 
thin specimens, there is not sufficient difference between the cathode-ray oscillo- 
graph traces obtained when the specimen is present and when the bars are wrung 
together for accurate values of the strain to be determined. On the other hand, 
if the specimens are too thick, the assumption that the pressures on both sides 
of the specimen are sensibly the same is no longer valid. 
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The method of determining the stress-strain curve from the photographic 
records is described in the section dealing with the theory of the method. The 
procedure after the plates had been measured was first to convert the measurements 
into displacement-time curves, using equation (20). The displacement-time 
curve obtained from the experiment when the specimen was present was then 
differentiated numerically, and this gave the pressure—time curve, as may be seen 
from equation (7). The strain-time curve was derived from the difference 
between the corresponding values of the two displacement—time curves, using 
equation (13). All measurements of time were made from the same point on 
the trace from the cylindrical microphone; as this microphone was in the same 
position along the bar in both experiments, this ensured that the values of stress 
and strain corresponded to the same instants of time. 

Before deriving the stress-strain curve from these results, the stresses had to be 
corrected for the radial kinetic energy effect, as described in an earlier section, 
and where the strains were large the change in the cross-sectional area of the 
specimen during straining had to be allowed for to obtain the true values of the 
stress. 

The experiments were all carried out at room temperature and calculation 
showed that the heating produced by the compression of the specimens never 
amounted to more than 1-5° c. in these experiments. 


(it) Polythene 

Figure 8 shows the stress-strain curves obtained with four different thicknesses 
of polythene (polymerized ethylene) of a type known as grade 20. The specimens 
were prepared from cast sheet which was carefully annealed to ensure that there - 
was no preferred molecular orientation. The points on each curve correspond 
to 2usec. intervals in the stress-strain cycle, so that the rates of stressing and 
straining may be derived from the curves. This same convention has been used 
in the dynamic stress-strain curves for other materials shown in the later figures. 
To avoid confusion around the origin in Figure 8 the first two or three points on 
each curve have been omitted. The broken curve on this figure shows the results 
obtained for the static stress- strain relation, using the apparatus shown in Figure 7. 

It may be seen that the effective elastic modulus for the dynamic tests is very 
much higher than that observed statically, and that in all four experiments the 
strain continues to increase after the stress has passed its maximum and has 
decreased to a fraction of its maximum value. Since there was no permanent 
plastic flow in the material, the specimens being of exactly the same dimensions 
after the experiment, this increase of strain with decreasing stress may appear 
at first sight rather surprising. The interpretation of this effect in terms of a 
memory function of the type first SueECaed by Boltzmann is discussed in detail 
in the next section. 


(iii) Rubber 

Experiments were carried out with specimens of natural and synthetic ribbed 
and the results obtained were somewhat similar to those observed with polythene. 
The delayed recovery effect with these materials was even more pronounced than 
for polythene, the strain continuing to increase rapidly after the stress had passed 
its maximum value. 

The specimens of the three synthetic rubbers investigated, neoprene GN, 
Perbunan and GR-S were closely similar in their static mechanical behaviour to 
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the specimen of natural rubber used, and it was found that their dynamic behaviour 


was also very similar to that of the natural rubber. 


For all these materials no 


‘permanent deformation could be observed after the dynamic tests. 
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Figure 10, Dynamic stress-strain curve for natural rubber. 


Figure 10 shows a typical dynamic stress-strain curve obtained with natural 
The specimen used for this experiment was initially 0-03 cm. in thickness, 


rubber. 
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and the maximum strain observed was about 20%. The broken curve shown 
on the same figure was one derived theoretically in terms of a memory function, 
and is discussed in the next section. No static stress-strain curve is shown, since 
on the scale used it would be indistinguishable from the strain axis, the static 
value of Young’s modulus for the material being about 2 x 10’ dynes/cm*. The 
difference between the static and dynamic behaviour of rubber is thus very much 
more marked than it is for polythene. 


(iv) Polymethyl-methacrylate 


This is a hard glass-like plastic of which Perspex is a plasticized form. Figure 
11 shows the experimental results given at 2sec. intervals for a stress-strain 
cycle, together with the static stress-strain curve obtained with a similar specimen 
in the hydraulic ram apparatus. The continuous curve in this case is the one 
derived theoretically and is discussed in the next section. 

It may be seen that the dynamic elastic modulus for this material is very much 
higher than that found with rubber and polythene, and that as the stress is removed 
the strain decreases, although the material does show considerable hysteresis 
at these rates of loading. 
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Figure 11. Dynamic stress-strain curve for Perspex. 


When thick specimens of Perspex were inserted between the bars it was found 
that a pressure pulse was communicated through them from the main bar to the 
extension bar. ‘The transmitted pulse was flatter in shape and of longer duration 
than the initial pulse, as well as being much smaller in amplitude. By using a 
series of specimens of up to 2 in, in thickness, it was found that the head of the pulse 
was propagated with a velocity of 2,400 metres/second; this corresponds to a value 
of Young’s modulus of 6-8 x 10!°dynes/cm?. In making this determination of 
the velocity, the separations between the head of the pulse and a fixed point on 
the trace from the cylindrical condenser microphone were measured on the 
cathode-ray oscillograph records. 

When comparable thicknesses of either rubber or polythene were tried no 
well-defined pulse was transmitted through them. This is to be expected from 
the dynamic stress-strain curves shown in Figures 8 and 10. From these it can 
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_ beseen that the layers of material will not recover immediately the stress is removed, 
_ so that a continued lengthening of the pulse with a corresponding drop in the 
| pressure amplitude will take place as the pulse travels through the material. 


(v) Copper 

Figure 12 shows the dynamic and static curves for specimens of pure annealed 
, copper sheet 0-05cm. in thickness. When the specimens were measured after 
the experiments it was found that they had decreased approximately 4°% in 
thickness. The behaviour of the metal is thus very different from that found 
/ with polythene and rubber, although the appearance of the curves is rather similar. 
_ The material is here flowing under the influence of the applied stress, the only 
) recoverable part of the deformation being the very small strain corresponding to 
| the initial steep rise in the stress-strain curve. The difference between the static 
} and dynamic behaviour for copper is not nearly so marked as it was found to be 
with the other materials investigated. 
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Figure 12. Dynamic stress-strain curve Strain 


for copper. ; ; : 
PP Figure 13. Dynamic stress-strain curve 


for lead. 


(vi) Lead 


Figure 13 shows the dynamic stress-strain curves obtained with lead specimens 
| 0-05 and 0-10cm. thick. As in the case of copper, it was found that after the tests 
| the specimens were thinner by the fraction indicated by the curves. As a result 
_ of the comparatively small stress required to start deforming the specimens and 
| the high density of lead, the radial kinetic energy correction was found to be large, 
and the values of the stress in the region near the origin were consequently rather 
uncertain. ‘The results did, however, show that the yield point was very small 
compared with the order of magnitude of the maximum stress applied in these 
experiments. 
| It was found that for this material the values of the stress were approximately 
| proportional to the rate of straining; this indicates that the lead is behaving like 
a viscous liquid under these dynamic conditions. 

As in the case of rubber the stresses required to deform the material statically 
were found to be too small to be shown on the same scale as the dynamic results. 
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§6. DISCUSSION OF RESULTS 


The experimental results described in the previous section have shown three 


types of mechanical behaviour that may occur when materials are submitted to } 


rapidly changing stresses. Under these conditions, polythene and rubber only 
begin to recover after the stress has fallen back to zero, but the recovery is complete 


and no permanent deformation can be observed after the experiments. Poly- . 
methyl—methacrylate shows a hysteresis loop at these high rates of loading, but : 
with this material a considerable fraction of the strain disappears within a few ° 
microseconds of the stress passing its maximum. Here again, after the experi- | 
ments, there is no permanent deformation and, when the specimens were examined | 


through crossed nicols, they showed no optical strain. 


Copper and lead, unlike the plastics and rubbers, were permanently deformed | 


and showed no recovery either immediately the stress began to fall or subsequently. 


The difference in behaviour between copper and lead results from the much higher ° 
yield point of the former metal, which has a definite elastic region in its stress— - 


strain curve. The yield point of lead on the other hand is so low that in these 
experiments its behaviour is very much like that of a viscous fluid. The static 
stress-strain curve for copper shows that even after the yield point has been passed 
the material continues to offer increasing resistance to further deformation as a 
result of work-hardening. 

In the interpretation of the results obtained with the plastic and rubber 
specimens it is necessary to assume that the relation between the stress and the 
strain is also a function of time. When a specimen of such a material is stressed, 
the strain does not reach its final value immediately but increases with time, 
Conversely, if it is deformed a fixed amount and the strain is then maintained at 
this value, the stress relaxes with time. 

Boltzmann (1876) was among the first to consider this type of behaviour 
analytically, and more recently several surveys (Slonimsky 1940, Leaderman 
1943, Alfrey and Doty 1945) of the subject have been published. Boltzmann 
introduced a new concept in his treatment of the problem, and suggested that the 
mechanical behaviour of a specimen was a function of its entire previous history. 
He further assumed that when a specimen had undergone a series of deformations, 
the effect of each deformation was independent of the others, and the behaviour 
of the specimen could be calculated by a simple addition of the effects that would 
occur when the deformations took place singly. This assumption has become 
known as the Principle of Superposition. 

In considering the behaviour of a specimen in which both the stress and the 
strain are varying simultaneously, we may consider the strains produced by known 
stresses, and the problem is then one of creep. Alternatively we may calculate the 
stresses produced when the specimen is deformed a known amount and treat the 
problem in terms of the relawation of stress. It has been shown (Gross 1947) that 
the two treatments are equivalent, and it is generally a matter of convenience 
which of them is followed. 

In Boltzmann’s original treatment he considered the problem in terms of the 
relaxation of stress and postulated that the stresses produced by dilatational and 
shear strains may relax in different ways. The case of uniaxial compression or 
tension involves both shear and dilatational deformations, and if two separate 
relaxation phenomena have to be taken into account the problem becomes rather 
complex. Boltzmann, however, suggested that it was unlikely that many materials 
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_ would show dilatational relaxation phenomena and in fact there is no experimental. 
_ evidence to show that such phenomena occur. Recent work (Pinkerton 1947, 


/ Mason et al. 1948) on liquids has shown that when such relaxation does occur,, 


) the times involved are very much shorter than those with which we are concerned 
| here. 

Boltzmann considered the past strain history of the specimen as divided into 
, a number of infinitesimal time elements. At time ¢ the residual stress due to a 
strain <(r) which took place at time 7 and was maintained for an element of time. 
| dv is expressed as — f(¢—7)e(7) where ¢(t—7) is called a memory function. If at 
i. time ¢ the value of the strain is <(¢) there will be a stress [«(z)] due to this, together 
) with the residual stresses due to all the previous deformations. Thus, the 
+ resultant value of the stress o(¢,¢) will be given by 


OBOE [ PEC An eS (24) 


In interpreting the delayed recovery effects observed in the present work on 
' plastic and rubber specimens, it has been found convenient to use a slightly 
1 different treatment from that suggested by Boltzmann. This treatment is 
) similar to that used by G.I. Taylor and E. Volterra (Taylor 1946) in their work. 
| onthe mechanical behaviour of polythene at lower rates of loading, and an analogous 
) approach has been described by Leaderman (1943) in his treatment of creep 
_ phenomena. 

| In this treatment, instead of considering the total strains to be maintained for 
| short elements of time, we consider the small changes in strain which have occurred 
7 inthe past history of the specimen and sum the residual elements of stress produced 
{ by them. ‘Thus, if between the times 7 and 7+ +67 the strain changes from e to 
| €+6c, this will result in a change of stress 5¢ which will relax withtime. If we call 
| the final value of the stress when relaxation is complete F(Se), then at time ¢ the- 
} residual value of the element of stress will be given by 


Bo=F(be)+f(t=7)8, wee (25) 


where /(t—7) is a memory function which approaches zero as ¢ increases. 
The resultant value of the stress o at time ¢ when the strain is ¢, is then given by 


o= Fl) + ( “ HOea) dr se (26) 


' A material whose mechanical behaviour is represented by (26) will not show any 
{| permanent strain but may show delayed recovery. 

If when the strains are maintained for a sufficiently long time the material 
& obeys Hooke’s law, F(e) will be given by Ee where F is the value of the appropriate. 
© elastic modulus. Taylor (1946) has indicated how the function f(t—7) can be 
determined from experimental results, but it was found that when this method 
was applied to the results of the present experiments the scatter was too large for 
+ asatisfactory determination of the function to be made. Instead, it was assumed 
) that f(t—r) was of the form A exp[—(¢—7)/«] and the values of A and « which 
» gave the best fit were then determined. 

\ This assumption is equivalent to regarding the mechanical behaviour of the 
« material as being similar to that of a system consisting of two springs in series, 
% with a dashpot connected across one of the springs. ‘The springs are assumed to 
obey Hooke’s law and the dashpot is assumed to obey Newton’s law of viscosity,. 
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i.e. the viscous force is proportional to the velocity with which it moves. F is 
then the effective elastic modulus of the two springs in series, and (£ + A) is the 
effective modulus for the spring which has no dashpot across it. The time of 
relaxation of the system is governed by the ratio of the magnitude of the viscous 
and elastic elements of the dashpot and spring in parallel. 

In deriving the theoretical curves which are shown in the figures the relation 
used has been 


d 
o-Be+A| exp[-(t-a)/al dr, asus 


the values of E, A and « having been chosen to give the best fit with the experi- 
mental results. It would appear possible to determine the value of E experi- 
mentally, since it should correspond to the value of Young’s modulus for strains 
which have taken place very slowly. Unfortunately, it is found (Alexandrov and 
Lazurkin 1940, Hillier and Kolsky 1949) that the value obtained for E with many 
of these plastic materials depends on whether the forces are applied in times of 
the order of minutes, seconds or milliseconds, and it is consequently difficult to 
determine the appropriate value experimentally. In theory, the value to be taken 
should be that for strains which have been maintained for an infinitely long time, 
but since experimental evidence shows that irrecoverable plastic flow as well as 
relaxation phenomena of very much longer period take place, the exact experi- 
mental significance of E in equation (27) is not so well defined, and it must be 
considered as an average value of the modulus for times long compared with the 
duration of the dynamic experiments but short compared with those of the static 
tests. 

The assumption that with these materials there is only a single relaxation time 
is a considerable over-simplification of the problem. ‘The assumption that 
Hooke’s Law is obeyed exactly with stresses which have been applied slowly is 
also not altogether justified. Nevertheless, fairly good agreement with the 
experimental results was obtained for polythene and polymethyl—methacrylate 
on these assumptions. With rubber the agreement was less satisfactory and it 
seems clear that the model assumed will give only qualitative agreement with the 
actual behaviour of rubber at these rates of loading. 

Figure 9 shows a plot of equation (27) with A =4-6 x 101°, H = 109, in dynes/cm2, 
and «=2ysec.; the values of « and de/dr were obtained from the experimental 
results with polythene and the observed values of stress are also shown in the 
figure. Using a similar treatment for results obtained with polythene, the stresses 
having been applied for times of the order of ten milliseconds, Taylor (1946) 
obtained lower values of A and EF and the value given for « is 1:7 milliseconds. 
Thus, it would appear that relaxation phenomena are exhibited over a range of 
times, but that in a given experiment the relaxation times of the order of the 
duration of the experiment are the only ones that are of importance. 

Figure 10 gives the theoretical curve for rubber, with 4 =2-4 x 10!°, F=108, 
and « is again taken as 2usec. The value cf F in this case is too small to affect 
the shape of the theoretical curve appreciably, and this value was chosen as one 
that corresponds to the modulus for stresses applied for times of the order of a 
few milliseconds. Although the values of A and « used were the ones which gave 
the best agreement with the experimental results, it can be seen that quantitatively 
the agreement is not very good and a more complicated model must be considered 
to account for the exact shape of the stress-strain curve. 
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Figure 11 gives the curve for d4=14:5x104, E=8 x10! and «=0-5 sec. 
‘The agreement with the experimental results obtained with polymethyl—meth- 
acrylate is here quite close and the difference in the behaviour of this material and 
that of polythene and rubber thus appears to be due to the shorter effective 
relaxationtime. This results in the term Ee becoming very much more important 
| in determining the instantaneous values of the stress than it was with either 
+ polythene or rubber. 


§7. CONCLUSION 


The experimental method described in this paper has enabled measurements 


} to be made of the stress—strain behaviour of materials in times of the order of 


= 20sec. The main disadvantage of the method lies in the fact that the specimens 
have to be compressed between two parallel surfaces, and that even when lubricants 
' are used, there is some uncertainty whether there is free lateral movement of the 
+ specimen, and whether the effects of friction between the specimen and the bar 
are not introducing large errors. Since, however, the use of different lubricants 
did not affect the results, and where permanent deformations were obtained they 
were fairly uniform with no evidence of barrelling, friction effects would not appear 
} to be large. In the experiments, it was attempted to maintain two thin lubricant 
_ layers between the surfaces of the specimen and the faces of the bars, and since 
) there was not sufficient time for the lubricant to be squeezed out, the conditions 
| are in this respect more favourable from the point of view of eliminating friction 
' than in comparable static experiments where the lubricant layers may become 
# very thin. 
| The interpretation of the results has shown that very marked relaxation 
_ phenomena occur in plastics and rubber in the microsecond region, and that their 
» mechanical properties under these rapid conditions of loading are very different 
{ from those observed when forces are applied more slowly. Although the 
_ assumption of a single relaxation time has led to reasonable agreement with the 
experimental results, the complicated molecular structure of these materials would 
| indicate that a large number of different relaxation phenomena take place simul- 
® taneously anda more complicated memory function is applicable. Work at lower 
# rates of loading (Alexandrov and Lazurkin 1940, Nolle 1947) has indicated that 
» the dynamic behaviour of materials composed of macromolecules is very sensitive 
} to temperature, and it would therefore seem desirable to continue these experi- 
| ments over a wide range of temperatures. 

The theoretical interpretation of rubber-like elasticity in terms of molecular 


| processes, which has been put forward by several workers (Guth and Mark 1934, 


Kuhn 1936), shows that the high extension and low elastic modulus of rubber- 
| like materials may be accounted for in terms of the uncurling of molecular chains 
under the influence of the applied stress; this involves the rotation of segments of 
the chains around single carbon-carbon bonds. ‘Two other molecular processes 
_ associated with very much higher elastic moduli may also occur in these materials ; 
these are changes in the values of the angles between the bonds and changes in 
| the distances between adjacent carbon atoms. Since the increase of interatomic 
} distances will require very large forces, it seems reasonable to correlate the 
/ two components of the elastic behaviour given in equation (27) with the changing 
) of bond angles and the uncurling of chains. Fis then the modulus corresponding 
|) to the latter process and the integral involving the memory function corresponds 
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to the fact that the uncurling process requires times comparable with the time of 
application of the stress. Thus, when the stress is first applied, the initial 
deformation corresponds to only a change in the bond angles. 
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ABSTRACT. The preliminary work on the measurement of velocity of propagation and 
attenuation of longitudinal sound oscillations in high polymer filaments already described 
by Hillier and Kolsky has been continued. It was concluded in that paper that temperature 
control would be advantageous, and that by varying the temperature and the frequency of 
oscillations more useful information could be acquired. An apparatus is described whereby 
the velocity of propagation and the attenuation of longitudinal oscillations throughout the 
range 0° to 50° c. and 500 c/s. to 30 kc/s. can be measured. The results obtained with one 
material, polythene, are discussed in terms of several theories of the elasticity of high 
polymers, and the constants of the several equations considered are calculated from the data. 
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long chain polymers when subjected to small, rapidly changing forces, in 

order to separate the elastic and plastic components of strain. The method 
employed was to send a sonic or ultrasonic longitudinal wave along a filament of 
the material, and to measure the velocity of propayation V and the attenuation 
per unit length «. The theory of the method used has been reported previously 
(Hillier and Kolsky 1949). The values of V and « were obtained for several 
frequencies between 500 c/s. and 30kc/s., and for temperatures between 0° and 
40°c. From the results obtained for polythene it has been possible to obtain 
values for an equivalent mechanical model in the manner described by Alfrey (1948), 
and to discuss other methods of expressing the results, (a) in terms of a complex 
compliance by analogy with the dielectric equations given by Fuoss and Kirkwood 
(1941) and others, and (4) in terms of Maxwell’s equation, derived from a definition 
of viscosity (Maxwell 1868). 


iE HE aim of the investigation was to study the mechanical behaviour of some 


S2) DESCRIPTION OF APPARATUS 


The apparatus can be conveniently considered in two parts: (a) the constant 
temperature enclosure containing the filament, and (d) the velocity and attenuation 
measuring equipment. | 

The constant temperature enclosure (Figure 1) consisted of a copper tube 
10 s.w.G., six inches in diameter and twelve feet long. This length was required 
so that measurements could be made on filaments whilst they were drawn to four 
or five times their original length. On to the outside of this tube was wound, and 
fixed by spot welding, a spiral of 16 s.w.G. wire with a pitch of three inches per 
turn. Over this was slid a second copper tube slightly shorter than the first, 
thus enclosing an annular space. ‘The wire diameter was such as to seal the gap 
between the two tubes so that when oil was pumped into one end of the annular 
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gap it was spiralled round the tube before emerging at the far end. The heat 
transfer was thus high, and by using a large, well-lagged reservoir tank and a high- 
speed oil pump the inside of the tube oven could be maintained at any constant 
temperature within the range. ‘T’o reduce heat loss and avoid local air pockets 
the copper tubes were mounted in a sealed wooden box of square cross section, 
so constructed that the tube rested in the centre with adequate air space all round. 
A centrifugal fan was coupled to the apparatus with flexible hose couplings, and 
air was driven along the outside of the tubes and returned to the fan via the inside 
space. The air system was closed so that both sides of the fan were used, and 
the pressure was nowhere very far from atmospheric. The oil pumped between 
the copper tubes was stored in a tank equipped with an immersion heater and 
refrigerating coil, and was lagged with 3 in. thick cork slabs. Suitable thermo- 
stats kept the oil temperature at a constant value in the tank, and the thermostats 
were sufficiently sensitive to maintain the air temperature in the oven constant 
to within 1°c. throughout the range 0° to 60°c. The temperature could be 
measured to 0-5°c. continually by means of a calibrated thermocouple. 


Copper Tubes 


Driving Motor 


Refrigerating Coil Pomp 


Figure 1. 


The measurement apparatus consisted of a crystal pick-up mounted on a 
carriage sliding on two long parallel cylindrical steel rods rigidly clamped together. 
This part of the apparatus was mounted in the constant temperature oven. The 
carriage was driven along the rods by means of a lead screw mounted between 
them and rotated from outside by a small reversing electric motor through a 
reduction gear and flexible drive. Provision was made for ascertaining the 
position of the carriage by using a flexible steel tape fixed at one end to the carriage 
and sliding past a fixed mark on the outside of the oven (Figure 2). The sound 
source was either a modified moving-coil mains energized speaker as described 
previously (Hillier and Kolsky 1949), used for frequencies up to 10 kc/s., or a 
magnetostriction transducer, used up to 30 kc/s. This transducer consisted of 
a nickel tube sealed at one end by a small pin chuck used to hold the filament 
coaxially, and with water-cooling tubes at the other end as shown in Figure 3. 


| 
| 
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‘The nickel tube was held in a ring clamp on which was mounted the energizing 
‘coil, Although no effort was made to drive the nickel tube at its resonant 
frequency, provision was made for altering the effective length to increase the 
response at selected frequencies. 

The crystal pick-up was connected to the main amplifier by a coaxial cable 
with insulated outer screen so that the pick-up, insulated from the oven, could be 
kept at the amplifier earth. The amplifier was a three-stage, wide band-pass 
unit with negative voltage feedback to stabilize it and was used with a Cossor 339 
double beam oscilloscope as previously. The sound source was fed from an 
R-C wide range oscillator (Dawe Instruments, Type 400B) through a suitable 
power amplifier. A valve voltmeter was used to measure the r.m.s. voltage on 
the Y, plates so that a measure of the amplitude variation along the filament 
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Figure 2. Pick-up carriage and drive. 
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Figure 3. Magnetostriction oscillator. 


‘could be made. ‘To obtain the velocity of sound from the measurements of 
phase shift the frequency of the oscillator must be known. A quartz crystal 
oscillator and frequency dividing unit were built so that 100 kc/s. and 25 kc/s. 
signals were available to calibrate the wide-band driving oscillator. ‘The quartz 
crystal used was of special low temperature coefficient cut, and it was considered 
that the frequency of 25 kc/s. was accurate to within 0:5%. As in the previous 
work, equipment was available to extend the filament at a constant rate, and to 
measure the tension in the filament at the same time as velocity measurements 
were being made. The apparatus for this is shown in Figure 6. 

For all the work at frequencies up to 10 kc/s. and for temperatures up to 40° c. 
an ordinary piezo-electric gramophone pick-up cartridge was employed together 
with the type of needle described previously (Hillier and Kolsky 1949). For 
higher temperatures and for frequencies above 10 kc/s. at all temperatures, a 
quartz crystal unit was designed so that the inertia effects, which reduced the 
output of the ordinary type of unit above that frequency, and the temperature 
effect, which rendered Rochelle salt useless above 45°c., could be eliminated. 

46-2 
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The unit is shown in Figure 4. For simplicity and ruggedness a simple bending 
system was employed; two rectangular silvered quartz plates 1 mm. thick and 
Sin. x Lin. in dimensions, cut with X axes perpendicular to the face, were cemented. 


Spring 


Filament 


Perspex 


Quartz 
Crystals 


le 
Figure 4. Quartz crystal pick-up. 


y Coaxial 
Cab 


together in a very light Perspex holder with a light copper spring used to grip the 

filament. The crystals were clamped at the end in a heavy brass clamp with 
thin rubber pads to protect them from damage. This unit was successfully 
employed for frequencies up to 30 kc/s. and with temperatures up to 60°C. 


§3. RESULTS 
The relation between observed phase shift and the propagation constants 
has previously been obtained (Hillier and Kolsky 1949), and if 6 is the value of 
phase shift at a distance / along the filament 


1 + me! 
tan 6 = ————1 an Rl 
1 —me-“! 
where m is the reflection coefficient at the pick-up, k is equal to 277/A, « the attenu- 
ation factor and A is the wavelength of sound oscillations in the material. The 


amplitude ratio B/A is 
B/A=(1—m)e—*/(1 — 2me- cos 2k1 + me— 8. 


For large values of J, 0=kl, and accurate values can be obtained from the slope 
of the graph of @ against /. Similarly for large values of J and «, B/A=Ke-™, and 
hence plotting log (B/A) against / resulted in a line whose slope gave the value of «. 
For intermediate values of « a different method was employed which was not so 
accurate (Hillier and Kolsky 1949). Another approximation is given by Ballou 
and Smith (1949). | : 

The method gave results of the same general accuracy as those described 
previously. ‘The frequency, however, was known to a greater accuracy and the 
temperature maintained more constant, so that more consistent results were | 
obtained. Experiments were carried out on a variety of materials throughout 
the range of temperature and frequency. The most comprehensive results were 
obtained on a sample of polythene supplied by I.C.I. Alkali Division, Alkathene 
grade 20. Annealed filaments approximately 0-025 in. diameter and 150 cm. 
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long were used. ‘Table 1 gives the values for velocity of propagation and the 
attenuation per unit length for this material. The results represent the mean of, 
several experiments having a variation of less than 3°% for the velocity measure- 
ments and 10°% for the attenuation. It is hoped to publish later further results 


_ with other materials and also the results for experiments when the material was 


drawn during measurements. ‘Static’ values, Ey, of Young’s modulus calculated 
from these stretching experiments with polythene are, however, included in 


(Lable:1. 
Table 1. Values of V (m/s.), « (cm™) and Eg (x 108 dynes/cm?) for Polythene. 


Frequency 


1 kc/s. 2 ke/s. 5 ke/s. 10 ke/s. 16 ke/s. 


—) Es V a V o Va os V a V a“ 


0 26:0 | 1041 0-005 | 1066 0-010 | 1108 0-019 | — = = SS 
10 10:92 | 866 0-005 | 903 0:013 | 926 0-032 | 944 0-047 | 960 0-086 
30 5730) 1) 608). 0-0084— 615 0-017.) 639..0-035. | 650. 0-077 747670 _ 0-114 
40 4°60.| 510-—0-007 |-520 0-022 | 535 -0-068 | 538 0-105 | 547 0-138 


§4. DISCUSSION OF RESULTS 


In the experiments carried out on polythene, the velocity of phase propagation 
V and the attenuation per unit length « for small amplitude Seu aie oscillations 
along the filament have been measured. 

In order to correlate these measurements with physical constants of the 
material, such as elasticity, it is necessary to consider the relation between stress, 
strain and time for the material, and then to derive the wave propagation constants 
V and «. 

The fundamental hypothesis put forward as a starting point in a study of 
the deformation of real materials is the superposition principle enunciated by 
Boltzmann (1876). A more general theory has been developed by Volterra 
(1912), which includes Boltzmann’s theory as a special case. ‘This theory, of 
considerable mathematical complexity, is based on the theory of functionals 
(Volterra 1931). ‘These are discussed in a recent paper by Kolsky (1949). 

From time to time equations have been proposed relating stress, strain and 
time for the deformation of. certain materials, which should have a limited use. 
Such equations are not intended to apply universally, but to be related to some 
specific conditions of material, stress or time (duration of stress). Maxwell (1868) 
suggested a differential equation to include both elastic and viscous behaviour, 
and since this equation has had very wide application and has formed the basis 
of much subsequent work, it is worth quoting. 

If o is the elastic force acting and « is the strain produced, then, for a material 
obeying Hooke’s law, o=e where E is the elastic modulus. If viscosity is 
present then o will vary with time at a rate dependent on o and the nature of the 


body. Consider this rate proportional to a, then 
do da. =o 
ee eet Sy Se 3 
dt E (| (3) 


1/7 being the constant of proportionality. 
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If de/dt is constant, o = Er(de/dt)+ Ce”. Then Eris the coefficient of viscosity 
and r is a time called the time of relaxation. A perfectly elastic body has an 
infinite relaxation time, and for a purely viscous fluid 7=0. Equation (3) does. 
not hold for many materials or over a wide range of stress and time intervals, but 
the idea of relaxation time has been extended and has proved very useful. 

The concept of separating the elastic and viscous components of strain has 
given rise to the use of models to represent the mechanical behaviour of materials 
and to allow for calculation of the results of specific applied stresses. The simplest 
such model is based on Maxwell’s treatment, where an elastic component and a 
viscous component are joined in series. Such models can consist of either 
mechanical or electrical components (Alfrey 1948, Gemant 1936). Either the 
behaviour of the material is calculated by analysing a system of springs (pure 
elastic elements) and dashpots (viscous elements) linked in some manner (mecha- 
nical analogue) or by computing the current or E.M.F. in an electrical circuit 
containing capacity and resistance (electrical analogue). Simha (1942) has given 
a peneralized model theory. 

Frenkel and Obrastov (1940) have published a generalization of Maxwell’s 
equation using tensor notation and the postulate of ‘volume viscosity’. Their 
final complete equation contains seven constants. The results of any ordinary 
tests made on real materials are very difficult to put into this very complicated 
equation, so that the evaluation of the constants and a verification of the equation 
are both difficult and tedious. 

Tobolsky, Powell and Eyring (1943) have published an account of a physical 
theory based on Eyring’s rate process conception and certain free energy 
hypotheses. In this theory a deformation is considered as affecting a kinematic 
unit of molecular dimensions, and for a macro-body the stress—strain—time 
equation is the sum of al] such small unit deformation processes. ‘Thermal 
agitation will result in an oscillation of the unit about its mean position in some 
form of ‘energy well’. If a potential gradient is formed by the application of 
stress, such that work is done when a unit moves in a specific direction, then the 
net velocity of units in the direction of the gradient can be evaluated by summation. 
When certain assumptions are made as to the kinetic units and the type of energy 
well in which they oscillate, this theory can be shown to be akin to the model 
theory. The treatment is of interest in that it incorporates temperature as a 
further variable and enables thermodynamic relationships to be developed. 

It is proposed now to consider three equations connecting stress, strain and 
time: first, Maxwell’s equation, secondly, one based on a model in which a spring of 
elasticity H and a dashpot of constant 7* are connected ia parallel, as in 
Figure 5(a), and thirdly, a composite model akin to the second, but with an added 
spring in series with the spring—dashpot unit:(Figure 5(5)). For each of these 
equations the constant of elasticity and relaxation time will be derived in terms 
of the velocity of phase propagation V and the attenuation «. 

In all cases w is the displacement along the wx axis, i.e. along the filament, and. 
a, the applied stress, is of the form o=oysinwt at x=0. 


(1) Maxwell’s equation. 


A de Co da du 
Fro t 3 aS eo ath PA ela 
m equation (3), Ee - + ape and=¢ ae 
Thus E d*u =e do 


The Measurement of Dynamic Elastic Constants 707 


By Newton’s laws, since acceleration is d°u/dt® and increment of force is do/dx, 
pd@?u/dt? =da/dx ig Oe HO) 


where p is the mass per unit length of filament. 
A solution of these two equations of the form 


U=Uyexpiw(t—x/c); o=o,expia(t—x/c) 
satisfies the condition c=o)sinwt at x=0. 


(a) | (8) 


Figure 5. 


—— Spring Balance 


Cord attached to filament 
Naw 


Cord tension shown thus a 


Figuce 6. 
The equation u =u, exp 7a(t — «/c) is a form of the wave equation with frequency 
w/27 and with complex propagation constant c=a+ib. Then 
U=Uy exp — bwx/(a* + b*) exp iw{t —ax/(a?+b?)} ...... (6) 


A damped wave propagated in the direction of x with frequency w/27 will 
give rise to a particle displacement u given by 


U = Uy Exp (— ax) exp iw(t—x/V). 
Thus the attenuation factor « =bw/(a?+b%) and the velocity V= (a? + b?)/a. 
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From equations (4), (5) and (6) 
a= \(E|p) {wr|(1 +0%r)} (1 +o%?)! +07} 


oe B= 4(E/p) for /(1 + e%r2)} {(1 + 222) — cos}. 
sae 7 =(w? — V02)/2V a0", } (7) 
E= V2ap|(w® sa Vx2), SOI ons 


Table 2 gives values of + and E corresponding to the values of V and « at the 
various temperatures and frequencies used. 


Table 2. Values of r (x 10~sec.) and E ( x 10! dynes/cm?) for Polythene 
calculated by Maxwell’s Equation 


Frequency 


Gee 1 ke/s. 2 ke/s. 5 ke/s. 10 ke/s. 16 ke/s. 

aie T E T E T 10; T E T E 

O- [0-950 4201 |2046 1-00" 80-24) toda pe 

10 | 1-15 0-69} 0-42 0-76 | 0-17 0-80 | 0-11 0-83 | 0-060 0-86 

30 | 1-0 0-3 | 0-47 0-35 | 0-22 0-37 | 0-099 0-39 | 0-065 0-44 
1-4 0-43 0-25 | 0-14 0-26 


(ii) Model, Figure 5(a). 

This model represents a system exhibiting a retarded elasticity and is different 
from the Maxwell treatmeat where two separate mechanisms of response to stress 
are postulated. Alfrey (1948) discusses the distinction between the two systems 
in some detail, but for the present argument it is sufficient to-note that a time factor 
corresponding to the relaxation time of Maxwell’s equation can be obtained for 
the model. ‘This retardation time is the ratio y*/E and can conveniently be 
represented by 7 as before. ‘The equation derived from the model, where o, as 
before, is the applied stress, < the strain produced, is 

= 4 de 
o=He+n ao. ) ee ee (8) 
and, as before, by Newton’s laws pd?u/dt2=do/dx. A solution of the form 
U = Uy EXp 1w(t — x/c) can be obtained where c=a+7b, and 


= 4(E/p) {(1 + w?7?)t =P 1) * b= 4(E|p) ‘(1 a5 w7?)t = 1}, 
whence = =-1=2Va/(w?-V0?);  E=pV?2w%(w? — V2a2)/(w? + V202)? 


‘Table 3 gives values of r and £ for this system. It can be seen from equations 
(7) and (9) that the product of the relaxation and the retardation time is equal 
to 1/w*, which furnishes a check on the computations. 
(11) Model, Figure 5 (6). 


The model shown in Figure 5(b) gives an explanation of the behaviour 
attributed to a retarded elasticity. It cannot explain the instantaneous elastic 
response or the fact that the velocity of sound at very high frequencies is still 
finite. In order to improve on this model a component to represent the residual 
elasticity, effective at high frequencies, has been included. To the system 
illustrated by Figure 5 (a) a second spring has been added in series. This leads 
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| tothe following equations where, as before, o and « are the stress and strain of the 
} total system and E, Ep and 7* are the constants of the springs and dashpots. 
_ Ife, and eg are the strains in each part separately so that « =«, +€,, then 


o = Eye, = Enen + 9*(dex/dt) 


_ and, as before, p(d?u/di?)=do/dx and t=7*/E,. The following differential 
) -equation connecting the displacement u and the constants of the system can be 


obtained: 
E,+Ep@u or du_ du du 
on Boedta. © Bat dvi dxdt. ¢ °° 2 _ 


This equation must reduce to Maxwell’s equation if E,=0 and to equations (8) and 
(9) if E,= oo. A solution to equation (10) is of the form u=wuy exp iw(t — x/c) 
where c is complex, and is given by 


ab yess P(E, + ER) + pERw*7? aos ipE,wr (11) 
cE ER(14+ 27?) Hdl oe) 
Now Os Pa 7 (12) 


Table 3. Values for 7’( x 10™3 sec.) and E’ ( x 10° dynes/cm?) for Polythene 
calculated from Model 5a. 


Frequency 


2 ke/s. 5 ke/s. 10 ke/s. 
al E , E / 


Equating the real and imaginary parts of 1/c? in (11) and (12), 


2a pwr 
a ee 13 
Va Ep(1+w*) ey) 
w'— Ve"  p(E, + Ex) +pEpw?7 


pee Var E,E,(1 +a?) 


For very long times, so-called ‘static’ loading, the response of the system will be 
given by the elastic response of the two springs in series, i.e. the ‘static’ modulus 
Eg is given by 1/E,=1/E,+1/E,. This modulus was measured in the course 
of the experiments carried out to investigate the change in velocity of sound on 
stretching the material, and values of this ‘static’ modulus at each temperature 
are given in Table 1. Then from equations (13) and (14), replacing E, and E, 
by Es, 

V w — Vo? 


T= > CO - 
2akg 2V aw 


Table 4 gives the values of 7 calculated by equation (15) at the various frequencies 
-and temperatures used. 
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As can be seen from Tables 2, 3 and 4, the values of F and 7 are not constant 
with respect to frequency, and this can be regarded as showing the existence of a | 
spread of relaxation or retardation times over the time scale. Experiments at — 
much lower frequencies carried out on high polymer materials of similar long chain * 
structure, e.g. rubber (Alexandrov and Lazurkin 1940), have demonstrated the 
existence of relaxation times of the order 0-1 sec. to 0-01 sec.; long term creep 
experiments (Leaderman 1943) have indicated that some relaxation processes may 
have characteristic times of the order of minutes or days, and very recently Kolsky 
(1949) has shown that at very high rates of loading, relaxation effects in times of 
the order of microseconds occur. Consideration of the molecular processes that 
can occur in a material ‘consisting of long chain molecules very slightly cross. 
linked or, as in the case of polythene, held by regions of crystallization (Bunn and. 
Alcock 1945), is consistent with such a spectrum of relaxation times. 


Table 4. Values of 7 (x 1073 sec.) for Polythene calculated from Model 5d. 


Temp. Frequency 
(°c) 
Dkkc/so 2 ke/se| oi ke/se |) 10 Ke/saate doikelss 
0 2:78 1-4 0:8 —- aa 
10 6:0 2°5 1-06 0-73 0-412 
30 5°5 2:6 1-4 0-626 0-439 
40 6:0 1:9 0-64 0-415 0-323 


It was mentioned previously that Tobolsky, Powell and Eyring had considered. 
a system based on molecular forces which introduced temperature as a variable. 
A brief summary of a simple system of this kind is given below, together with its. 
application, to give an approximate estimate of the energy of activation of the 
process responsible for the retarded elasticity effect in polythene at the frequencies 
and temperatures given in Table 3. k 

The molecular unit is assumed to exist in a potential well, of form 
F=3AF(1— cos 2nx/A), where A is the distance between minima, x is the distance- 
from the unit and AF is the free energy of activation, that is the height of the 
barrier. The effect of an imposed stress field will be to add a linear term —c-Ax 
where a is the stress and A the eftective area of the unit. 

The total free energy as a function of x is F=4$AF(1— cos 2zx/A)—o Ax. 
Tobolsky, Powell and Eyring have shown that when cAA<2kRT and oAA> F,. 
then 

de MAY do: MAG 


A 5QARE Ge Re 


This can be identified with = = oe + 50 
1 wA me AF 
h -—- = = = _- 
aie E \2mAF? »* Ah RT 
a : da A a : 
Of, it 77 /E = 7, Wie aos (cf. equation (3)). 


eee ee 
Then T= 5 = 5a qpexP AF IRI, ie., Inc =K+AF/RT 
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Hence, if one molecular process is involved in the retarded elasticity effect, 
| Inz plotted against 1/T should give a straight line of slope AF/R. 

From the results can be obtained the order of the activation energy for the 
| process involved in this high rate low amplitude deformation. The value obtained 
| for activation energy is approximately 200 cal/gr. mol. at 2 kc/s. 

It has been suggested that a different method of expressing the results might 
# lead to a more useful theoretical treatment. In the case of a dielectric exhibiting 
_a frequency variation of dielectric constant and power factor the results of experi-- 


03 


0 


At Re peal : 
4000 6000 8000 10000 
Frequency (c/s) 


{000 2000 20000 


! Figure 7. Values of the ratio of the imaginary to the real parts of the complex compliance against- 
log frequency. 


1947). 


» Kirkwood 1941). 


/ the compliance. 


U =U, Exp w(t —x/Cc), 


where c=a+ib, it can be shown that 
pC’ =(w®— V2x2)/ Vw; 


pe! 


Table 5. Values of C’ ( x 107° cm?/dyne) and C” ( x 10744 cm?/dyne)* 
Frequency 
Temp. 

(°c). 1 ke/s 2 kc/s 5 ke/s. 10 ke/s. 16 ke/s. 
La GC (OL (Ce Co Cos (Gis (EM (Ox Gt 

0 0-988 1-65 0-941 1-61 0-874 1:18 —- —|]— — 
10 1-44 ~— 1-99 {232-2 -48 Nie?) ee Som ples let eile OB, 
30 2:95 4:60 2°88 4-8 2°68 3°84 | 2-58 4:14 | 2-43 3-73 
40 4:24 4-84 4-06 7-46 8°82) 5: 9075)| 3250" 6°88!) 5°05 55°50 


* C’ and C” are the real and imaginary parts of the complex compliance for polythene. 


| ments can be considered in terms of a complex dielectric constant (Fuoss and 
y By analogy it is possible to consider a complex compliance 
(reciprocal of elastic modulus) on a similar basis (M.I.T. Acoustics Laboratory 


In a simple material pc? =1/C, where c is the propagation constant and C is 
Then if C=C’ +7iC”, we have p(C’ +7C’’) = 1/c?. 
Now if particle displacement u is given by 


aoa. 
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where, as before, V, « and w are velocity of propagation, attenuation factor and 
pulsatance respectively. In the case of dielectrics, the ratio of imaginary to real 


{ 


parts of the complex constant is equal to tan6, the loss in the dielectric. Then . 


C’/C’ can be regarded as representing the mechanical loss in the material. 
‘Table 5 gives values for two quantities C’ and C”, and Figure 7 shows the variation 
of C’’/C’ with frequency for one temperature. 


§5. CONCLUSIONS 


It is apparent that the description of the mechanical behaviour of polymers — 


and the prediction of such behaviour in specific circumstances of load, load rate 
and temperature present very great difficulties to the physicist. ‘The several 
methods discussed in this paper are not exhaustive and several more have been 
discussed in the literature (Alfrey 1948, Slonimsky 1940, Houwink 1937, 
Leaderman 1943, Simha 1942, Tuckett 1942, Astbury 1949). 

Those mentioned were selected as most likely to apply to the experiments 
‘carried out. The values of the elastic constants obtained will hold only for 
‘conditions of stress loading similar to those met with in the experiments, 1.e., 
small amplitude forces applied for times of the order of milliseconds; both 
vibration damping and acoustic transmission are obvious practical applications 
of such forces. From the calculated free energy values it appears that quite 
small units of the molecules are involved in the transmission of these forces, 


chain segments of the order of 5 to 10 carbon atoms only being involved. This | 
conclusion is, however, very suspect, since the results indicate that the supposition _ 
of a single operative relaxation process is unjustified. Any experiment can only» 
select relaxation effects that have constants in its immediate time scale, so that 


altering the temperature will merely select a different set of relaxation processes 
for study. To regard the results as implying that the processes involved at 
temperature A respond differently at temperature B is inaccurate. Until more 
details are known about the spread of relaxation effects, calculations of activation 
energies for the processes involved will give very low values. 


The general conclusion must be made that an elucidation of the problems | 
of high polymer molecular physics requires that very wide ranges of both time 


and temperature should be investigated before a solution can be reached. 
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ABSTRACT. 'Thetheory of boundary waves (Young, Rubinowicz, and others) is developed’ 
for the two-dimensional case and applied to discuss image formation for small apertures, 

obstacles and phase-retarding laminae. The theoretical predictions are illustrated by~ 
diffraction photographs. It is hoped thus to contribute towards a fuller understanding of 

microscope images for biight-field, dark-field and phase-contrast conditions. 


§1. INTRODUCTION 


N 1917, A. Rubinowicz, starting with the well-known Kirchhoff scalar- 
| diffraction theory, showed that the effects of diffraction by an aperture in a 

thin screen could be written down as the sum of (a) the amplitude which could 
be predicted if the domain of the direct light were strictly limited by rectilinear 
propagation from the source, plus (b) the integral of certain amplitude contri- 
butions derived from the boundary of the aperture. 

This theory can be applied simply to explain the edge-diffraction effects which 
occur in important types of image formation, especially dark-ground illumination 
and phase-contrast. ‘Though quantitative calculations are convenient only in 
simplified cases, the idea of the ‘boundary waves”’ gives a more comprehensible 
account of the phenomena than can be obtained from the usual discussions given 
by Rayleigh and others in terms of sine and cosine integrals. 

The equation given by Rubinowicz for the amplitude w is as follows: 


1 il exp (tkp,) exp(tkr,) cos (n, 7.) sini(pndends 


=Uy =U agra 
tesla tle Uy m7 Ps (ee ECOs(7, (p,) 
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where uy is the amplitude of the “direct”, and up that of the diffracted light; 


the symbols in the integral expression may be understood by reference to Figure 1 


‘Figure 1. Figure 2. 


-where L is the source, F is the diffracting aperture and B the boundary; P is the — 


point at which the amplitude is taken. The element ds of the boundary is at a 
distance p, from L and 7, from P, 1.e. 


— od 

Eds=p.. Pde. 
and, finally, 2 is the inward directed normal (at ds) with reference to the cone 
defined by rays from L which touch the elements of the boundary. 


§2. SEMI-INFINITE SCREEN 
Referring to Figure 2, the direction of the incident parallel light is the Z axis 
of Cartesian coordinates, while the opaque screen occupies the X, Y plane below 
the X axis; thus the direction m is represented by OY. We consider a point P 


in the Y, Z plane defined by the polar coordinates 7, ¢ (positive as shown by the 
arrows). ‘Then 


(ei) Sr shecs a yo 


and cos(r,, p,) = — Roos f 
ie i. 


The Rubinowicz integral therefore becomes in our case 


: 
vaigt ie [_ SED AGP) 7 SID Pac pas ait ak re (2) 
re ees Yr,  1rcosd—-r, 
where A is written for the amplitude of the plane wave in the XY plane. 
Approximate numerical evaluation of this integral has been carried out. 
It leads to results similar to those obtained by the usual discussion of the Fresnel 
integrals and the associated phase-vector diagram (Cornu spiral), except in the 
immediate neighbourhood of the shadow, where complex considerations are 


required. We shall therefore now confine the discussion to the two-dimensional | 


case. It will be seen below that the expression of the resultant amplitude as the 
sum of “direct illumination amplitude” plus ‘“‘boundary wave amplitude” 
characteristic of the Rubinowicz expression is applicable in the more familiar 


treatment. This view was explained in a previous paper (Martin 1943), and a more 
detailed numerical discussion will now be given. 


{ 
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§3. THE DARK-GROUND IMAGE OF A STRAIGHT-EDGE 
The amplitude and phase in a plane surface S at a distance / from the diffracting 
edge E (Figure 3) and perpendicular to the incident parallel light is represented 
by the length and direction of a line such as G,F, (Figure 4) joining a point in 
the spiral curve to the focus. The coordinates C and S of points in the spiral 
are given by 
C=|cos(m0*/2)dv and S= [sin(aro?/2) do, (3) 
where Me a Aly et hw he Fhe ens Bat Montesa (4) 


and y is the distance of the point (represented by the value of v) within the 
- geometrical shadow. The coordinates of the focus are C=S =}. 


S 


Figure 3. Figure 4. 


Now since the gap between the circular locus C centred in E (Figure 3) and 
the plane surface S is (for small values of y) given by y?/2J, the associated phase 


difference is wy?/Al, that is, 7v?/2. Hence the phase ¢ in the circle C at the point 
having the same value of y is 


@=arctan {($—S)/($-—C)}-—av7/2, ss ss se (5 a) 
and the amplitude A is given by 


A? =(4—S)2+(4-C)2. 


Since the angle of slope of the spiral curve is 7v?/2, the value of ¢ is the angle 
between the element of the curve and the line joining that element to the focus. 
Starting from a value of 7/4 when v=0, it will asymptotically approach 7/2 as 
v increases. ‘Thus, subtracting a constant amount, 7/2, from these values, we 
may exhibit the relative phase angles as starting with a value of — 7/4 and approach- 
ing a constant value of zero. ‘Table 1 is calculated with the aid of the known 
Fresnel integrals. 

Apart from the region where v is small, the change of phase becomes slow; 
and we may regard the disturbances as approximating to waves (the boundary 
waves) centred near the diffracting edge. Ifthen the disturbance passes into an 
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image-forming system of small aperture, the above amplitude and phase distri-- 
bution can be regarded as transferred to a corresponding locus, APR (Figure 5),. 
centred in the geometrical image point E’ in the image plane.* Ap and ¢p are 
the amplitude and phase respectively at P, and «’ is the angle AE’P. Then the 
components of the amplitude at Q due to an element of the wave front at P are 
proportional (within well-understood limits of accuracy) to 


Apcos {¢p+(2mh'/A)sina’} and Apsin {6p +(27h’/A) sin a’}, 


where h’ =E’Q, being very small. | However, these expressions can be simplified 
for numerical computation. If R is the extreme point of the transmitted wave 


Table 1 
; Relative : Relative 
yp — Ampiide phase (fr) o Amplitude phase (de) 
(radians) J (radians) 
0-0 0-707 —7/4 225 0-140 —0-060 
0:25 0-560 —0:571 2°50 0-126 —0-048 
0:50 0-437 —0-410 2-715 0-115 —0-042 
0:75 0-352 —0:297 3-0 0-106 —0-035 
1-0 0-287 —0-218 2-25 0-098 —0-030 
1-25 0-243 —0-167 3°50 0-092 —0-025 
1-50 0-208 —0-122 S275 0-087 —0-023 
1°75 0-178 —0-092 4-0 0-082 —0-020 
2:0 0-157 —0:075 | 


Figure 5. 


front, and the angle AE’R is «/,, let us write (assuming radian units for W,) 
We= (Cah simon ee (6) 
then W,= (2zh' /A) sing’ “Wa (Upite), «5 5S ace (7) 


where vp and vp are the v values corresponding to the elements of the entrant 
wave front associated with Pand R. Note that W, will change in sign for points. 
Q above and below the axis. 

The intensity J due to the wave AR is obtained from the equation 


l= Gen ee eee (8) 
where S =k| A peeks + W(vp/vp)} dv 9 
S oP sin tPF R\Up/ UR en Me PR So cma (9) 


k being a constant. The phase of the resultant is arctan S/C. 


: The assumptions are of an obvious character, and need not be discussed in a treatment of 
the present kind. 
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Numerical case. 


In order to illustrate the procedure, the range of v is taken from 0 to 4, and 
numerical integrations are carried out with the extension of Simpson’s rule which 
is expressed by 


2 
J ‘ me (79 + 32y, + 12y + 3293 +7y4), 


where «x is the common interval between the ordinates yo, y1,....y. The effect 
is to find the area under a curve of the fourth degree laid through five points on 
the true curve which have equal increments x of their abscissae. This rule is, 
however, applied four times in the whole interval. If the ordinates were all 
unity the total area would be given as 16. It is therefore convenient to take x 
as 10/16, so that the full sum with unit ordinates should be 10. 


Table 2. Relative Displacements due to Lower Part of Boundary Wave only 


W C: S C24 S? W Cc S C257 
14 0-384 0-216 0-194 = il 1-809 = ilezHe 4-885 
13 O23 72 0-178 0-170 == 7) 1-268 —1-666 4-384 
12 0-353 0-173 0-155 = 0-676 seule Seal 
11 0:326 0-268 0-178 = A Ons7 = eats 2345 
10 0-393 0-349 0-276 ==5 5) 0-040 = 1 ODNS 1-498 
9 0-512 0-371 0-400 = © 0-022 = 1015 1-030 
8 0-614 0-316 0-477 = 7 0-052 =O3e 0:844 
7 0-632 0-237 0-455 = 0-021 943 0-889 
6 0-584 02375 105397 = 9 —0-083 = O81 0-873 
5 0-575 0-361 0-460 —10 —0-196 = Oso 0-755 
4 0-726 0555 0-835 = ili =S29s) ==e7Al 0-570 
3 1-092 0-662 1-630 Sa? 02229 = (Pave 0-404 
2 1-578 0-511 AST SN =iS — 02174: — (0-542 0-324 
1 1-992 0-046 3-970 14 —0°145 —(0:-541 0-292 
0 POUT —0-626 4-706 


Ke 
W, (RADIANS) ——> 


Figure 6. 


The result of the computation is given in Table 2 and Figure 6, in which the 
intensity distribution characteristic of the edge-diffraction is represented. It is 
the halo, or dark-ground image, which would be obtained if a screen (suggested 
by S,, Figure 7) could intercept the direct light while allowing all the diffracted 
light to enter the objective. There are important differences between this result 
and that which would be found (suggested by the broken curve, Figure 6) if the 
edge were self-luminous. The maximum intensity does not coincide with the 
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geometrical image position, and the darker bands which border the central 
maximum are not well developed. However, if the image-forming system is 
arranged so as to admit no light except that diffracted at reasonably large angles, 
the dark-ground “halo” image will approximate more closely to the type of 
the ‘‘self-luminous ” one. 


§4. DIFFRACTION DUE TO BOTH PARTS OF THE BOUNDARY WAVE 


If, instead of producing the dark-ground effect by the screen S,, in front of 
the objective, we remove the direct light by a small screen, S, (Figure 7), in the 
back principal focus, the boundary wave which spreads from the diffracting 
edge into the directly illuminated region can also take part in the formation of the 
final image. Referring to the vector diagram of Figure 4, the total effect at a 
point in the region above the geometrical shadow is typified by G,F;, which 
can be regarded as the resultant of G,F, and F,F,; that is, of the boundary wave 
represented by G,F, and the direct light F,F,. This is identical with the 
Rubinowicz principle. This new boundary wave G,F, has the same values of 
the phase for the same numerical values of v* but has the opposite sign for the 


Figure 7. 


amplitude. In the previous section we were dealing with the boundary wave 


—> —> 
G,F,; we now deal with G,F,, and it is obvious that when v is given any particular 
numerical value (for example when v = 0, and both G, and G, coincide at the origin), 
these amplitudes will have opposite signs. This is in accord with the modified 
Rubinowicz expression (2). 

Thus the components of the disturbance at the point Q due to the sum of the 
boundary waves without the direct light component are 


C UR COS oR cos 
5 =k i Ap ae {dp + Wp} du +k |, (— Ap) ee {bp — Wp} dv, 
that is 
oR OR 
C= —2k iF Apsingpsin Wpdv,... S=2k [ Ap cos gp sin Wy dv. 
~ 0 


Each of these integrals will thus have a numerical value which is equal but of 
opposite sign for points in the focal plane at equal distances above and below 
the axis. ‘The relative resultant amplitudes, i.e., (C2 +.S2), are shown in Table 3, 
and in curve 1 of Figure 8; the relevant phase angles w appear in curve 2, and the 


* vy represents the numerical distance of a point in the Cornu spiral, measured from the origin 
along the curve ; thus symmetrically placed points have the same v values. Y 
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intensities in curve 3. ‘The image of the edge should show a perfectly black line 
bordered by bright fringes if the image-forming light is limited to that of the 

boundary waves; this is easily observed in practice if a very narrow stop is placed 
in the principal focus of the lens in a system similar to that specified above. Plate I, 
A, shows the normal image of a straight-edge, and also the dark-ground image. 
‘The latter shows the centre dark fringe, bordered by bright fringes, as predicted 
above. Details of the system for the experiment are given below. 


Table 3 
Relative Displacements due to both Parts of Boundary Wave only; 
v=0 to + 4:0 

WwW C WO C+S27(C?24+S%)F ow W (a S C?+$2(C?+S8%)t w 
d(xadians) (deg.) | (radians) ( deg.) 

16:0 0:503 0:782 0:865 0-930 57-2 6° 0:562 1-252 1:883 1:372 65-8 

15:0 0-515 0:805 0:913 0-956 57-4 540) 0:535 1:584 2-795 UNIX, giles 

14:0 0-529 0-757 0-853 0-924 55-1 4:0 0:489 2:068 4:516 2°125 76:6 

13:0 0-546 0-720 0:817 0:904 52-7 3-0 0-415 2-383 5-851 2-419 80:25 
B62 61270) «690-582 0°766 0:925 . 0-962 Sy 2:0 0:310 2:177 4:835 2:199 81-9 
> 11-0 0-579 0-979 1-294 1-137 59-4 1:0 0-183 1-316 1-765 1-329. 82:8 

10:0 0-589 1-196 1-777 1-333 63°8 0:5 0:082 0:687 0:479 0:692 83:2 

9:0 0:595 1:302 2-049 1-432 65-4 0:25 0:042 0:352 0:126 0:354 83-25 

8:0 0-593 1-259 1-937 1-392 64°83 0:0 0-0 0-0 0:0 0-0 — 


UE) SoU Jol Om i719 eles 63-7 


~ -_--- 


AMPLITUDE 
PHASE 2 
INTENSITY 3 


EFFECT OF BOTH PARTS 
OF BOUNDARY WAVE 


Figure 8. 


§5. IMAGE OF A SMALL APERTURE 


Consider a small aperture G,G, (Figure 9), axially symmetrical with regard 
to the image-forming system, and illuminated by a parallel beam of coherent 
light ; the amplitude in the plane surface containing the pupil of the system is now 
given in terms of the Cornu spiral by the vector, G,G, say (Figure 4), representing 
the total interval v corresponding to the aperture. When we consider the central 
point of the pupil, the points G, and G, are symmetrical with regard to the origin, 
but they travel along the curve (maintaining a constant interval) when we estimate 


47-2 
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the amplitude at points further and further from the centre. We can distinguish 
three cases: 
(a) G, and G, both in upper part of the spiral: then put G,G,=G,F, + F,G2; 
(b) the origin lies between G, and G,: then put GG, =G,F, + FF, + F,G5"% 
(c) G,and G, both in lower part of spiral: then put G;G,=G,F, + FG. 
The components of the light entering an aperture of the objective system _ 
may therefore be described thus: 
(i) both parts of the boundary wave corresponding to G,; 
(ii) both parts of the boundary wave corresponding to G, (opposite sign) ; 
(iii) the direct-light component entering an aperture of the objective equal to 
the aperture in the diaphragm. 


Figure 9. 


4- | AMPLITUDE 


Figure 10. 


It is not difficult in some cases to calculate the light distribution in the final 
image. We have already calculated the amplitude and phase distribution in the 
image corresponding to both parts of the boundary wave; one such distribution 
is NOW centred in the image of G, while another, with opposite sign, is centred in 
the image of G, as suggested in Figure 10. We must further calculate the dis- 
placement in the image plane due to the direct light and then combine all three 
components with due regard to phase. 

In order to deal with any particular case arelation must be established between 
the above scale of distance (W,) in the image plane and distances h perpendicular 
to the axis in the object plane. 

Assuming that the pupils of the objective coincide with the principal planes 
and that h,h’ are the object and image heights, while the conjugate den 


ee Cees a both direct-light and boundary-wave contributions to the amplitude 
g ric shadow-edge, noted by Rubinowicz (1917), are manifest in this formulation. 
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are J and l’, we have h/l=h'/Il' provided that the system has object and image in 
air. Now Wp =(27/A)h'y,,/l’, where y,, is the semidiameter of the pupil; then, 
since, from (4), J=2y?,/Av2,, we find 


Lt eo), teers el Beek, MY 2) isthe (11) 


There is the complication that the distribution of phase calculated for both 
parts of the boundary wave does not apply immediately to the normal image plane ; 
this will be obvious from Figure 9. The calculated phase applies to an image 
surface inclined at an angle tan-4(h+h’)/l’. The relevant correction must be 
made. We may, however, ignore the correction to the numerical value of the 
amplitude and also ignore for a first calculation the fact that the boundary wave 
calculations made above for an axial ‘“‘edge”’ should really be corrected for an 
object position off the axis. Both these corrections are likely to be very small 
in our case. 

With reference to the direct light, its phase in the entrance pupil is uniform ; 
therefore its advance over the phase of a wave centred in the axial point of the 
object plane is wv?/2. On transmission through the optical system this wave 
(of uniform amplitude) will maintain the same phase-lead over a convergent 
wave centred in the axial point of the image. It is therefore simple to carry out 
integrations for any numerical case provided the phase changes concerned are 
not large. 

These principles can now be applied in a numerical case, and the one chosen 
should lend itself to experimental observation. 

Let 7= 1,000 mm., /’ = 200 mm. and v,,=4. TakeA=0-5 x 10-3? mm.; then the 
corresponding value of y,, is 2-0 mm. (equation (4)). Let the breadth of the object 
be 0:25mm., each edge being at a distance of 0-125mm. from the axis; the 
corresponding image distance is 0:025mm. We now have the weg corresponding 
to the image distance from the equation (11); itis radians. Moreover the value 
of v corresponding to the edge of the direct-light beam in the pupil is 0-25. 

‘The calculation for the boundary-wave image will apply to a plane which 
intersects the conjugate point and makes the angle of (0-125 +0-025)/200 radians 
with the usual image plane. The corresponding phase corrections may thus be 
computed ; their signs will be clear from Figure 9, but it should be noted that the 
distance /’ is greatly Paoeseated in the diagram, in comparison with the actual 
case of calculation. 

With regard to the direct light (amplitude=A, and phase dyp, say), using 
the same notation as in equation (10) above, but adding a suffix d to signify the 
direct light, 


cos ¢, 
2 veh 2kA ain fa 608 (Wa 2) ell (12) 


and in accordance with previous provisions 
dpa SS av/2 — m/4. 


Aq is taken as 4/2, the value of F, F, in the vector diagram. 

The relative displacements due to the “direct light” are given in ‘Table 4. 
The results are corrected so that they have the proper ratio to the figures of 
Tables 2 and 3. 
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The results of the numerical calculations are given in Figure 11; the “dark- 
ground” curve shows the results of the boundary waves (Figure 10) added to- 
gether with the relevant phase corrections, while the curve marked “aperture”’ 
shows the result of adding the direct light; the broken vertical lines at w= +7 
show the limits of the geometrical image, and it will be obvious that this image, 
of which the breadth is twice the resolving limit of the objective for self-luminous 
object points, has no geometrical similarity with the object. 


Table 4. Relative Displacements* due to Direct Light (v=0 to +0-25) 


W Amplitude W Amplitude W Amplitude 
0 ie y 3-0 ileys 8-0 1-68 
0:25 1°77 4-0 1-74 9-0 1-66 
0-5 1:76 5:0 1-74 10-0 1-64 
1-0 1:76 6:0 1°73 11-0 1-62 
2:0 1:76 7:0 ile7l 1270) 1-60 


* The variation of the resultant phase angle (—47°) is negligible for present purposes. 
Note that diffraction has spread the illumination beyond the geometrical image. 


ae 
10 APERTURE 
/ 1% DARK-GROUND 


Figure 11. 


§6. SMALL OBSTACLE 


‘The case of image-formation when the object is a small obstacle obstructing | 


the light is dealt with in a very similar manner. ‘The amplitude in the pupil 
is represented by the interval F,F, between the end points of the Cornu spiral 
less the interval G,G,, which represents the interval of v corresponding to the 
obstruction. In the same three cases as before we have 


(a) F\F, —G,G,=F,F, —(G,F, + F.G,); 


(d) =F\F, -(G,F, + F,F,+F,G,); 
(c) = FF, - (G,F, + F,G,). 
These components now represent the direct light amplitude from the whole 
aperture minus the boundary-wave amplitudes of the previous case. We have 


therefore to add the direct light effect for the full aperture (with the appropriate 
phases taken into consideration) to the original amplitude contributions taken 
with reversed signs. According to Lummer’s principle, which can be invoked 
when the aperture is assumed to be unrestricted, the direct light will provide a 
background image of coherent light of uniform phase; this assumption is not 
strictly valid in our case, but it may be adopted as a first rough expedient. We 
will therefore assume that the amplitude (of the background light) is constant 
and the phase angle is zero (but subject to the above conventional correction of 
—7/2). As for the numerical amplitude value of the direct light, it will be of the 
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same order of magnitude as the values given in Table 4, though not so weakened 
by diffraction. The exact magnitude need not be calculated for this rough 
discussion, and we will therefore assume that a likely value is about 4-0. 

The remaining curve in Figure 11 shows the bright-field image intensity 
distributions thus calculated. The edge-diffraction light will cause a series of 
weak fringes which border the dark image of the obstacle. The dark-field image 
of the obstacle is identical with that of the narrow aperture, as far as intensity 
is concerned. 


§7. IMAGE OF PHASE-CONTRAST OBJECT 


The effect of passing a plane wave through a phase-changing lamina may be 
described as composed of ‘the original wave plus the wave due to the object”’. 
In Figure 12, if the original wave (amplitude and phase) is represented by OA, and 
the effect after transmission by the lamina is OB, then the amplitude and phase 
of the ‘‘ wave due to the object” is AB. 


eal 
oO A 
Figure 12. 
Table 5 
Relative Amplitudes* and Phases in the Image of the Small Aperture 
+We 0 1 2 3 4 ) 6 7 8 9 10 


Amplitude Deo ap-O4) S20 7°30 GedlsO2 "0-148 NO-62 260-45. 0:14" (0825. > 0°22 
Phase (dea) 06; 97 97) 89" 88" —90 — —-69° 70 —73 T t 


* Intensities are plotted in Figure 11. 
+ Phase calculation here subject to large uncertainty owing to small values of 2C and 
uS. The ratio &S/2ZC gives the tangent of the phase angle. 


The problem of the image formation of such an object is thus related at once 
to the problem of the aperture, and this has already been treated in terms of the 
‘“boundary waves plus the direct light”’. 

To illustrate the application, let us suppose that the object is a thin uniform 
lamina of the same size as the aperture and obstacle above, and further, that the 
retardation is such that AB/OA is 1/5 and that AB then shows retardation of 
(roughly) 7/2 with respect to OA. We can therefore take one-fifth of the values 
of the “‘aperture” amplitudes (Table 5) to represent AB, and the phase angle can 
be put as roughly 90° throughout, while the relative amplitude and phase of the 
direct light (assumed uniform) can be taken as (say) 4-0 and zero respectively. 
Figure 13 shows the resulting curves for ordinary conditions, and also for phase 
changes of +7/2 for the direct light, to secure positive and negative phase contrast 
(lower and upper curves respectively). Though the visibility of the “‘image”’ 
of the phase strip is greatly enhanced by the use of phase contrast there is, of course, 
no particular gain in the rendering of its geometrical form. ‘The case calculated 
by Francon (1946) is that of a wider object; the edge-haloes which he predicts 
are comprehensible in terms of the boundary waves. 
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§8. EDGE HALOES 


We havealready noted that cases arise, especially in the case of the microscope, 
where the object plane contains many small obstacles having sharp boundaries ; 
some may be more or less opaque, and some may be edges of, or steps in, phase- 
changing laminae. These sharp edges will diffract their pairs of boundary waves 
inio the objective and there will be, in the symmetrical case, a corresponding 
amplitude distribution of the type of curve 1, Figure 8, in the image plane, in 
addition to the effect of the direct light, which may be considered to give uniform 
amplitude. But the relative phase of the direct light will depend on the inte- 
gration of all the contributions derived from the object-plane, and for various 
reasons it will not usually be possible to predict it. A variety of bordering fringes 
may therefore appear according to the relative phases of the boundary wave system 
and direct light contributions, and scme of these are suggested in Figure 14. 
This argument applies only to cases in which the angular aperture of the 
illumination is very small. 


Figure 14. Intensities of edge-halo fringes according to relative phase of background light. 


A. Retardation 90°. B. Retardation 45°. C. Retardation 0°. 


§9. DIFFRACTION PHOTOGRAPHS 


The photographs B (Plate I) and C (Plate IJ) illustrate the general conclusions 
for the dark- and light-ground images of apertures and obstacles (slits and wires). 
The illumination was from a long-focus collimator. Green light (A=0-5 
approximately) was used, the width of the slit diaphragm in front of the objective 
was 4) mm., the distance of the object 1,100 mm., and the distance of the image 
ZS WM. 
ahe Images so produced were projected by an eyepiece magnifying x 2-49 
and the resulting negatives were enlarged with a further magnification of x 5-70. 
Uhe arrangements are thus fairly near to those of the case for which calculations 
were made above. In the case of the opaque wire object of diameter 0-262 mm., 
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the illumination in the central part of the image corresponding with the curve of 
Figure 11 is clearly seen. The image of the aperture of similar size shows no 
sign of any darkening at the centre, and the bordering fringes are very faint 
indeed compared with the border fringes in the dark-ground image. The interval 


| corresponding to the “maximum to minimum” of the (sin W)/W distribution 


corresponds to 0-148 mm. in the obstacle. The dark-ground images were pro- 
duced by stopping off the greater part of the direct light by a very narrow wire 
placed in the principal focus of the objective. It is seen that (as expected) the 
dark-ground images of wires and slits are indistinguishable. When the size of 
the object, be it slit or obstacle, is changed, the effects in the dark-ground image 
can easily be predicted by a corresponding change in the separation of the points 
G,,G, in the amplitude curves (Figure 10). 

Starting with a broad object, the points Gj, G3 are widely separated, and wesee 
two intensity systems of the type of curve 3 of Figure 8, one corresponding to each 
edge. As the separation diminishes, however, the nearer bright bands coalesce, 
and it will be found that the amplitude distribution approximates more and more 
closely to the type (sin W)/W, the position of the first dark minimum at each side 
of the central maximum depending only on the aperture of the system; however, 
the amplitude will fall continuously till the points G and G, coincide, when there 
is, of course, complete extinction both in theory and practice. It is a very inter- 
esting thing that the amplitude curves of the type of curve 1 of Figure 8 should 
have this particular additive property, though of course the “(sin W)/W” type 
of distribution is predictable from the Cornu diagram. 


§10. CONCLUSION 


The foregoing discussion has been given in the hope that it may afford a starting 
point in a more suggestive and useful approach to the interpretation of the micro- 
‘scope images of small objects in cases where a severely restricted aperture for the 
illumination is employed to obtain dark-ground or phase-contrast effects. It has, 
unfortunately, been considered by some writers that Zernike’s (1942) phase 
vector theory of phase contrast can be applied directly to elicit information about 
the relative optical thicknesses in the close neighbourhood of loci of phase-change, 
-or even in phase-modulating objects of sizes not much exceeding the ordinary 
resolving limits. Such information is not likely to be obtained from a static 
photograph, but it does not appear impossible that a good deal of information 
about the object can be obtained by controlled variation of the amplitude and 
phase of the “direct light”? and study of the corresponding diffraction effects at 
the boundary. ‘The discussions in this paper afford no more than a beginning; 
the results cannot be applied directly, for example, to images secured by means 
of phase-contrast with ring illuminators and the like, but similar principles can 
be used, and it is hoped in future to extend the calculations to a greater variety of 
-cases. 
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ABSTRACT. The wave-front aberrations which occur in a pencil traversing a symmetrical 
optical system ata finite field angle are defined in a new way. Refraction and transfer 
formulae are derived for all the aberration coefficients up to and including those depending 
on the fourth power of the aperture of the pencil. 


Sh INTRODUCTION 


N a symmetrical optical system the wave-front aberration, or deviation of 
the wave front from an arbitrarily chosen reference surface, is a function of 
three variables p,, ¢,, o defined as follows: 


Wave Front 


Principal Ray 


Figure 2. 


Let A’ (Figure 1) be the centre of the exit pupil, A’O’ the paraxial principal ray 
to the Gaussian object point O’, and P’ a point in the exit pupil. Then p,= AUP 
o is the angle between A’O’ and the axis, and 4, is the azimuth of A’P’ with 
respect to A’O’, 1.e. the angle between the meridian plane and the plane containing 
P’ and tA'O". 

It may be shown (Hamilton 1833) by considerations of symmetry, that the 
wave-front aberration W’ at P’ can be expressed as a power series in which the 
variables are powers and products of p2,p,0cos¢,,02. When this is done the 
terms fall into groups of the fourth, sixth, eighth, . ... degree in p, and o, these 
being the primary or Seidel aberrations, and secondary, tertiary, .. _, aberrations. 
The Seidel aberrations are those terms containing p‘, p2a cos ¢,, p20? cos 26 
par cos ¢,o%, i.e. of the fourth degree in p, and o. sub fond io 

However, if all the terms of degree 4 or less in p, are considered, irrespective of 
their degree in o, a number of new terms are brought in, containing p?cos°d 
Pa cos “has p,,cos*¢,, multiplied by various higher powers of o. This micas that 
if the field angle of a pencil is large, aberrations of these types will be appreciable, 


and since they are only of the third and fourth degree in p, they will be important. 


Auth 
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even in pencils of low aperture. Since these are not Seidel or primary aberrations, 
to obtain their coefficients in the usual way would involve computing secondary, 
tertiary, and even higher order aberration terms, depending on how important 
the higher powers of o were in the system concerned. Dr. H. H. Hopkins, to 
whom the preceding argument is due, has suggested that a better approach would. 
be to suppose that a true principal ray has been traced trigonometrically through 
the system, with s- and t-fans, at the desired field angle, and to determine the 
aberration coefficients up to the fourth power of the aperture of the pencil 
surrounding this ray. These coefficients would be expressed in terms of the 
principal ray data, and they would give the wave-front aberration in the form 


W’ = Boop” cos 2h + 
B;,p? cos@ + Bs3p? cos 3d + 
Biop- + Byop*cos*p+ Bap*tcos*@ = = wiae.. (1) 


referred to the sagittal focus. In equation (1) W’ is the optical path between the: 
wave front and a spherical surface centred on the sagittal focus and touching the 
wave front where it meets the principal ray. W’ is taken as positive when the 
wave front is in advance of the sphere; this corresponds to the usual sign conven-. 
tion for longitudinal ray aberration. Byo, Bs,,.... By, are aberration coefficients. 
p and ¢ are the aperture radius and azimuth in a plane perpendicular to the 
principal ray and meeting it at some suitably chosen point, e.g. the intersection 
with a refracting surface or with the axis of the system. 

Of the terms in (1) those containing B,5, Bz, and By are of the same kind as. 
Seidel astigmatism, coma, and spherical aberration respectively but not numerically 
equal to them, because of the finite field angle, while the remainder are of higher 
order types, although still only of third or fourth degree in the aperture. 

This procedure has the advantage that all variations in the aberrations with the 
field angle are taken into account, since the principal ray is accurately traced. 
Thus equation (1) is accurate up to the fourth power of the aperture, no matter 
how great the field angle may have been. In the usual procedure for studying the 
aberrations of off-axis pencils at finite field angles, the principal ray is traced, 
with s- and ¢-fans, then a pair of meridian rays at finite distance from the principal 
ray and possibly also a true sagittal (skew) ray. However, it can be seen that the 
data thus obtained will not differentiate between coma of the type p?cos¢@ and 
that of type p*cos*¢ since, from (1), the wave-front aberrations due to these 
two vary in the same way in the meridian and sagittal sections (6=0 and 6=7/2), 
and the difference appears only in intermediate azimuths. Similarly no differen- 
tiation between types p*cos*¢é snd p*cos*¢ is possible. On the other hand, if 
these coefficients could be computed directly, a complete picture of the behaviour 
of the wave front in all azimuths would be obtained, accurate to the fourth power 
of the aperture. This was pointed out by Dr. Hopkins. Formulae for the 
aberration coefficients will be derived in this paper. 


§2. DEFINITION OF THE ABERRATION COEFFICIENTS OF A PENCIL 


The term “ pencil”’ is used to denote the ensemble consisting of a true principal 
ray (not a “‘ paraxial”’ ray) from an object point at a finite or infinite distance to- 
gether with all the wave fronts, or surfaces of constant optical path, and all the 
rays, or wave-front normals, from this point. The aperture is, of course, to be 


728 W. Weinstein 


thought of as limited in extent, but the exact degree of limitation is not material 
to the discussion. One can then refer without further explanation to the wave 
front at any point of the principal ray or to the ray from any point on a certain 
wave front, and s0 on. 

{t is convenient to compute aberration coefficients which are simply related 
to those defined by equation (1). 

Take right-handed rectangular axes with the origin at some point of reference 
A on the principal ray, chosen as in §1, the x-axis along the principal ray, and the 
y-axis in the meridian plane (Figure 2). If the principal ray has traversed the 
system at finite angles of incidence the wave front at any point will have the meridian 
plane (xy plane) as a plane of symmetry and in general it will have no other kind 
of symmetry. Hence the equation of the wave front at A may be written 


1 


1 nee 1 
x= rea i 9,2 + Cao” + Crays* +r (a + Cu) 9 


1 1 
oF (a +Cu) ye (53 + Cu) ce 
+ terms of degree 5 and higher. 


In this equation only terms with even powers of z have been included, to give 
the required symmetry about the x—-y plane. ¢ and s are the radii of curvature 
of the wave front in the meridian and sagittal sections, at the origin. In §§4, 5 
and 6 below, the origin is taken at the intersection of the principal ray with the 
refracting surface, and then ¢ and s are equal to the intersection lengths determined 
by the ordinary formulae for tracing ¢- and s-fans. Cy, Cyp,.... Cog are the 
aberration coefficients which will be determined. The terms in s and ¢ which 
are included in the coefficients of y*, y2z? and z* ensure that, when all the C’s 
vanish and t=s the equation will represent a sphere, to the fourth degree in y 
and z. Itis not essential to include these terms in ¢ and s, since (2) is the Cartesian ~ 
equation of the wave front, and however the aberration coefficients are defined 
they must, when substituted in (2), give the actual shape of the wave front; 
nevertheless, it simplifies the algebra to define them so that there is no aberration 
when they vanish. Whenever higher order aberrations are defined the definitions 
are variable to a certain extent. Thus in (2) the coefficient of y2z? could be written 
as 1/(4ts”) +Cyp instead of 1/(42?s)+Cy, but, as explained above, if Cj, were 
computed instead of C,, the actual shape of the wave front would still be known, 
so that no uncertainty would arise. Such a change from C,, to C3, would merely 
mean that a slightly different definition of primary astigmatism was being tacitly 
adopted, the difference being in the shape the wave front is assumed to have 
when only primary astigmatism is present and the fourth power of the aperture 
is taken into consideration. 

To obtain the actual wave-front aberration from (2), referred to any chosen 
focus, it is simply necessary to subtract from (2) the equation of a sphere centred 


on that focus and to multiply the result by the refractive index N. (Thus a 
sphere centred on the sagittal focus would be 


I 1 
Ho = 7 (yi +22) + so(ta7yPt.... 
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and the wave-front aberration referred to the sagittal focus would be N(x —x), 
where N is the refractive index.) This gives the aberration as a function of y and 
%, rectangular coordinates in the exit pupil. If the substitution y=pcos¢;. 
2=psin¢ is made in (2) the right-hand side becomes 


delay ot u 
5) (7 ~ :) p* cos® $ + 3,P + (Cs — Cy2)p® cos? ¢ + Cy2p® cos ¢ 
iLjAll 2 1 
siP {Cio — Con + Cor =(3 = Rs aia 3) } pteost¢ 


iy Sa | ‘| 1 
=F {C22 ae ila = 3) \ pt cost b+ (Co.+ 5a) p* 


On subtracting the corresponding expression for x9, multiplying by N, and. 
comparing this with (1) it can be seen that 


‘ed aa 
Bu=>(5- 3): Bs =NCjp; 


Bgx = N{ C39 — Cio} 5 By =NCou; 
Accents 3 
Buy=N{Cu~2Co.t 74 - 3) | ©) 


4|t?s 3° 


(By ae | 
Bu=N{Co-Cor+ Cou (a “Wigs 02 a) | 


Peis, 
These equations relate the coefficients defined in (2) to the more usual polar 
coefficients in (1). 

It can be seen that the aberration as defined by N(x —x ) is an optical path 
length measured parallel to the principal ray. Aberration is usually thought of 
as measured either along the ray through the point P or along the normal to the- 
reference sphere, but in the former case the direction would depend on the amount 
of aberration present and in the latter on the choice of reference sphere. The 
present definition is independent of these factors. The difference is small with 
ordinary numerical apertures, but the point must be settled for the definition of 
higher order aberrations, and for this purpose the choice given here is most 
convenient. 


§3. GENERAL METHOD OF DERIVATION OF THE ABERRATION 
COEFFICIENTS 

Methods of deriving aberration coefficients generally depend on the principle 
of obtaining a quantity which is invariant on refraction, expressing this as a series 
of powers of some measure of the aperture, the coefficients in the series containing 
the aberration coefficients, and writing down the equations which state that 
each term in the series is invariant on refraction. ‘These equations are the re- 
fraction formulae for the aberration coefficients. Finally, transfer formulae, for 
transferring the coefficients from one surface to the next, are found; these are- 
usually fairly simple. The measure of the aperture must be a refraction invariant, 
and it is usually the central angle or some closely related quantity (von Rohr 1920). 

The method used in this paper follows similar lines, and may be described. 
as follows. 
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Let AP (Figure 3) be a refracting surface, AO the incident principal ray, 
AO’ the refracted principal ray and AQ the incident wave front at A. The 
refracted wave front AQ’ can be found in the following way. Let Q, Q’ be chosen 
so that the normals at these points pass through P and let N, N’ be the refractive 
indices of the two media. Then Huygens’ construction gives N! Q’P=N. QP, 
ort O P=¢g OP =9¢, 

ACNG) =O.) 0 See? eae cee (4) 


Part of the refracted wave front at A so constructed is in most cases “ virtual” 
(as in Figure 3), but this is immaterial since the real wave fronts at points further 
along the refracted principal ray are all derived from it; the aberration coefficients 
of this wave front will be found, and the refraction invariant * which forms the 
basis of the method is the quantity Ng. 

If A is taken as origin, the coordinates of Q will be (x, y, x) where x is given in 
‘terms of y and z by equation (2), since Q is on the wave front whose equation is 
(2). The coordinates of P will be, say (u, v, w), where uw, v, w are functions of 7, 


Figure 3. 


‘the radius of the surface, the distance AP, and %, the angle between the meridian 
plane and the plane containing P and the normal to the refracting surfacef. 
A numerical value of # is fixed and inserted in the expressions for u, v and w; 
-in what follows this will be 0, 7/4 or 7/2. 
Let (J, m, n) be the direction cosines of the ray through Q; these can be found 
as functions of y and z from (2), by analytical geometry. Then the fact that the 
ray through Q passes through P is expressed by the equations 


m=arhlggy sy) © 9 ae aera (5 a) 
Vay Te nige Be 77 ea eee (5 d) 
WSN | Be ee ee ee ee (5 c) 


where x, /, mand nare expressed in terms of y and z, and u, v, w are in terms of AP. 


_* The reason for using equation (4) is that it is valid for all positions of P, not merely when 
P is in the meridian plane, and hence it can be used to obtain aberrations in any azimuth. Suppose 
- for the moment that P lies in the meridian plane, and let PQ meet the principal ray in a point O, ; 
then 4N(AO,-PO,) is the aberration introduced at the surface, referred to the point QO, as focus, 
and this expression can be used to obtain all orders of tangential aberrations ; but, since there 
is no intersection point O, when P does not lie in the meridian plane, this expression is useless 
for other aberrations. Again, if O. is the focus for close sagittal rays 1N(AO,-PO,) can be used 
to derive primary aberration expressions (Hopkins unpublished), but at finite incidence angles it 
is not permissible to assume that skew rays pass through the s-focus to the required order of accuracy 
even for aberrations varying as the third and fourth powers of the aperture. 
+ The angle 7 is not equal to ¢ defined in § 1, unless there is normal incidence. 


a 
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From these three equations y and z can be eliminated and an expression for 
-q obtained ; this will be a power series in AP, and the coefficients will be functions 
of the principal ray data and the aberration coefficients Cy9,.... Coy. This 
expression for g and a similar one for q’ are substituted in (4) and the terms in 
the successive powers of AP are separately equated to zero, giving explicit ex- 
pressions for the quantities A(N cos IC,,), A(N cos 9JC39), etc. These are the 
refraction formulae (equations (13), (28), (14), (29), (20)). 

The incident wave front is a “ parallel surface ”’ to the refracted wave front at 
the previous refracting surface, i.e., it is obtained by laying off a constant distance 
D along the normal to the refracted wave front. This can be done by analytical 
geometry, and the result is a set of formulae giving the coefficients of the new 
wave front (aberration coefficients) in terms of those of the old. ‘These are the 
transfer formulae (equations (36) to (40)). 


§4. REFRACTION FORMULAE WITH #%=0 


This case is the simplest, since the rays all lie in the meridian plane, and the 
problem is only two dimensional. 


P Refracting Surface 


Principal Ray 


It 


Figure 4. 


Let AP (Figure 4) be the meridian section of the refracting surface, AO the 
incident principal ray, AC the normal to the surface and AQ the meridian section 
of the incident wavefront. Let QPO be the ray of the incident pencil which passes 
through Q. Let, be the radius of the refracting surface, J the angle of incidence 
of the principal ray, and ¢ the distance from A to the tangential focus for rays 
close to the principal ray, so that ¢ is the radius of curvature at A of the curve AQ. 
Take rectangular axes with origin at A and x axis along AO. Then, from (2), 
the equation of the curve AQ is 


Seis xy + C9 V2 + (a - Cw) y*-+terms of higher order, ......(6) 


the terms in z being omitted since, as stated above, only rays in the meridian plane 
are considered in this case. 
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Let AP be denoted by p; then the coordinates (wu, v) of P are expressed in 
terms of p as follows: 

Take for the moment x’, y’ axes with origin at A and x’ axis along AC, the 
normal to the circle AP. The equation of this circle is then x? +’? =2rx’ =p?. 


Hence 


, 1 / 1 
x paae VED rae 


MePS ao’ +O(B”), 

where O(p”) denotes terms of the mth and higher degrees in p, the aperture of the 
pencil. ‘The value of g which is finally obtained as described in §3 must be correct 
to the fourth power in #, i.e. terms in g which are O(p?) can be ignored, since 
aberrations depending on the fourth power of the aperture are being investigated. 
By following through the algebra the orders of magnitude of the terms in g to 
which any intermediate expression contributes can be seen, and any non- 
significant terms in this expression can be dropped. ‘This omission of non- 
significant terms will usually be carried out without further explanation in the 
remainder of the paper. 

Changing back to the x, y coordinate system, 


u=psin I+ cos VS 2 sin I+O(p°*), | 


a Pease ciel 
v=pcosl 7, Sin I Acie 


These are the required coordinates of P, expressed in powers of p. 

To find the direction cosines (/, m) of the ray QP, equation (6) is considered 
to be written in the form F(x, y) =0 by transferring all the terms to the left-hand 
side. Then 


ENG}: ENG By 


OF Cis, 171 
From (6) xe me Ls ae Feo 3 CoV? + O(y) ; 
OF \2 OF \2 1 6 
1 S: 
Then =1- apy 7 C209” + O(y4), | 
; Fo eae eee (9) 
Sagres OW), 


These are the required values of the direction cosines, which, with the values 
of u and v given in (7), are substituted in (5 a) and (55) giving 


Pr : 1 1 
p sin I+ > cos I— Fin = 598 + Cus?+ (5 + Cw) 


1 Spe 
+441 — irs 7 Cy" + Ol D9). iene eee (10 a) 


‘Dae : 1 
pcos [— 7, Sin I- Ecos Tay 91F9+3Cp9} OUD tae terse (10d) 
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It is next necessary to eliminate y from these equations and solve the result 
for g. From (104), by the usual method of reversion of series, 


3 
y= =p cos I~ F sin 14 cos 1- Ecos I+ (30s cost I~ 5 sin 1) pg 


22 EAC em a Tha ee (11) 


This is substituted in (10a) and the sci equation is solved for g, again by 
reversion of series, giving 


Posin re4 ( ee co) 
Gi 2 
+ {¢sin nt (Se - oe) - gasin I C55 cos® ihe 


ri cosI cos? 2 ay is cosf oc cos? I lee te 
+ {a ; pei. ) + 5g sin i Z =e 372, Sin® i, 


Perc 1(222 s <=) Cap — Cyy 6084 1} p+ 0(0%, 


This value of g is substituted in (4): since p is invariant on refraction the 
various powers of p can be taken outside the A sign, and since (4) is true for all 
values of p the coefficients of the various powers of p in (4) must each be zero. 
This gives 

A(N sin I) =0, Snell’s law. 


cosf cos? I 


av( ) =0, the equation for tracing f-fans. 


Ty 


NsinI /cosI cos?I \ 
o7\ _ 
A(NC gq cos? J) af ay ( : ? ) ; 


2 . 2 2 
A(NC 4 cos* 1) -an {3 sin I cos 1(< Hi ) Chet: (= ee! ) 
i % 


t r t 
drew COS, = 4cOsed. ities 
rr 1 ae ) = gasin ds eco (14) 
Equations (13) and (14) are the required refraction formulae for Cy) and Cy. 


A result equivalent to (13) was obtained in a different way by Dr. H. H. Hopkins 
and was communicated privately to the author. 


§5. REFRACTION FORMULAE WITH #=7/2 


The usual coordinate system is taken with the x axis along the principal ray. 
- The point P is taken on the great circle of the refracting surface which is tangent 
to the z axis at the origin A, i.e. in azimuth 7/2 with respect to the normal to the 
refracting surface; this great circle lies in a plane inclined at an angle J to the 
zx plane (see Figure 5), and it can be seen from the figure or from the analysis. 
below that neither P nor Q lies in the zx plane to the order of approximation 
contemplated. The coordinates of P are 


tenes Rt ip = 2 Of 15 
u=7 cost; v=—Z int; w=p—gatO(p’). (15) 
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The coordinates of Q are (x, y, 2), where x is given in terms of y and z by (2), 
but since y is of the same order of magnitude as 2, 1.e. y is O(p?), it is permissible 
to put 


1 1 1 
w= 5 B+ Cypya® t+ ay" + (x3 + Cu) A+ O(p?). sees (16) 
As in §4 
uman (Se, > Se)/{(Ge) + (5) + (G5) 
sm n= (E Oy’ Oz / Ox oy Oz ry 
so that 
ie ee 4 

l= 1 = le He O(p 15 
m= OPN 5m pooh 2 ee ices (17) 


Sr 


n=—+s+O(p%). 


7 Principal Ray 


~ Refracting Surface 
Figure 5. 


Substituting from (15), (16) and (17) in (5a), (5) and (5c), 
2 


p 3, ea ae 1 1 
77 C8 L= 5,8 + Cyyz er re seu (53 + Cu) 2+4(1- 3524) + O(p°), 


“a 


Pa cs (18 a) 
Bane I=y+O(p4) (183 
Dr ya ‘PJ, ee act eterersts ) 

p® 1 . 
pga a= ee OG ies ee (18¢) 


In these equations q is taken as O(p”); this follows from (19) below. 
Substituting for y and z from (185) and (18 c) in (18a) and solving for g by 
reversion of series, 


a 7 . ih COST ING ‘lata a 
+( (3 9S ea eer ee + gasint I(5 — 2) 


Cre. 
ry >, sin Ue Gn. pt ODS) ie eee ee (19) 


l| 


q 
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Applying (4), the successive terms give 


AN (= — *) =0, the equation for tracing s-fans, and 
A(NCo4) = AN] S#sin 14 La (= = } + asin? 1(- re :)}. 


‘Equation (20) is the refraction formula for C,,; it involves the values of Cj, 
before and after refraction, as might be expected, and the refraction formula for 
this coefficient will be obtained in the next section. 


§6. REFRACTION FORMULAE WITH ¢=7/4 


In this case P is taken on the great circle through A, the plane of which is in 
azimuth 7/4 with respect to the norma! to the refracting surface (Figure 6). 


f 


RA 


Refracting 
Surface 


Figure 6. 


The coordinates of P can be found as in § 4, by means of the auxiliary axes 
Ax’, y’, 2’. For typographical convenience put a=1/,/2. Then the x’, y’, 2’ 
coordinates of P are 


ods oF 3 og po 
Gypit as ByOES ar + O(p°). 


Hence in the x, y, 2 system, 


_ ng ap? Z 
u=ap sin I+ 7, or (p*), | 
Drs ap cos 1~¥ sin I~ F cos 14+ O(P%), Cee (21) 


———EE— 


wW=ap— =: p?+0(0°). 


736 W. Weinstein 


The equation of the wave front must now be taken as in equation (2), since all the — 
terms given there are of significant orders of magnitude. The direction cosines 
of the ray QP are, retaining all significant terms: 


1 3 1 ‘ 
et a 7? : Czy? — (7 sie 2) Cuays? +0104), | 
a= ry — SCagy" = anee as O(p*), | eisistsiers (22) 
1 | 
n= — 5 2 2Creya + O[f"). | 


Substituting from (2), (21) and (22) in (5a), (54), (5c), 


2 ; ae i ‘tae 
ap sin 1+ 4 cos I cE sin [= oe 3,2 + Cao 


i 1 1 
+ Cyoy2? + (3 + Co) 9+ (a +x) yrs + € +Cos) gt 


1 1 3 ey 
+441 Sor ae Go (; + =) Cray + O(P”), 


scene (23 a) 
pe. ap? 1 ey. 
ap cos I—"7 sin I- 3,3 008 I=y-q FY + 3 Coy” + Crp3" + O(p'), 
ES in re ae (23 b) 
ap- giPi=2—q{=242Cyyel + O14), re (23 c) 


As before, y and z must be eliminated from these equations, and the result 
solved for gin terms of p. It can be seen from the diagram or from (25) and (26) 
below that ; 


y=ap cos 1+ O(p?), g=ap sin I+ O(p?), 


so that in the highest degree terms which are significant in any equation it is. 
permissible to put y=ap cos J, g=apsinJ. Using these, it is found that 


{a a 


a : 
z= apt —pgt 138 sin? [ — 3,5 +aCy, sin I cos ihe OCD) a ge dees (24), 


2 


y=ap cos 1—Fsin I+ = pq cos [+ {jpesin® J cos I- <a cos I 


a a 


tet . 
~ Dy_ Sin I+ 5 Cyq sin I cos? J+ 5 Cy, sin Th ps + OD?) Tee. (25) 
J: 
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The values of y and z from (24) and (25) are then substituted in (23a), 
| neglecting terms which are O(p°), and the result is solved for g, giving 


cos?J °° 1 


q= ansin r+ {Ot — SE hye Foon Fin I~ sin 7 
<5 cos? ee eee ae TLS Cys cos Th 


+ 2 ine 3 ae in2 
ee cost [ — SB cos*/ sin? sd: Srp cos? [ + cos J sin? I 


4rt? 
+ eR coe? I- fa 3, im [+ ae : xasint I- se cos J + He 
—— 2 cos? J+ ? sin J cos (et ~ — t) C30 
cece {sin I (- — ee a =) Cyp — 4 Cg cos* 1 
—4C,, cos? J — Cou} Pe Os) An ee EO ea eis (26) 


Using (4), the first term on the-right hand side of (26) gives Snell’s law, and 


the second is a linear combination of the formulae for tracing s- and t-foci. The 
third term gives 


AN{ Cyp cos! 14 Cyg 08 1+ nt BE ae? 2Tsin I— ~ cos Isin I =); 


2s" 27. 
wares (22) 
If equation (13) is subtracted from this, the result is 
A(Ci2N cos 1) = AN sin I cos J — =asin if. Snel . (28) 


Equation (28) is the refraction formula for the aberration coefficient Cy, 


_ The term in p* in (26) gives in a similar way, on using (24) and (20) to eliminate 
i, and-C,,. 


D : 
A(C,,N cos? I) = AN} oe (+ _ 2cosI 408 ‘) C, 


t 
ifs sin I a ease si sin? f= ! at ee te 
Dr 0 ou sie 2rt? me oh 4rt? 4724 
: jee 2 fie - iE cos? [ (29) 
ANE abe 5,3 81m TG 73 COS T Ay Reoneiseete 


This is the refraction formula for Cys, the last to be found. 


$7. TRANSFER FORMULAE FOR: THE ABERRATION COEFFICIENTS 


All that now remains to be done is to find transfer formulae giving the values 
of the aberration coefficients before the pencil is refracted at a surface in terms of 
their values immediately after refraction at the previous surface, as explained 


in § 3. 
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Suppose that the equation of the wave front immediately after refraction at a 
certain surface is 


1 , 1 ’ Peat ND Sef 1 , / ; 
= hy + fet+Chy/?+ Cay's'*+ (goa +Ci) 9” ; 


1 1 / / 
cts (aay ais Cin) y't9"*+ (a +i) Ae As ere eee (30) 


and let the equation of the wave front immediately before refraction at the next 
surface be (2), i.e. the same as (30), but with the primes omitted. Let the distance: 
along the principal ray from the first surface to the second be D’; then if points 
are taken at a constant distance D’ from the surface (30) along its normals, the. 
locus of such points will be the surface (2) (Figure 7). 


x’ 


Principal Ray 


Refracti 
Refracting Surface Surface S 


Figure 7. 


If the normal to (30) has the direction cosines (/’, m’, 2’) the points concerned 
will have coordinates x, y, z given by 
x=x'+(l'—-1)D’, 
y=y +m'D’, aos Ce bea (31) 
g=2'-tn'D’, | 
But since the two surfaces have a common normal, /', m’ and n’ in (31) may be 
replaced by J, m, n, the direction cosines of the normal to (2), giving 


x’ =x—(1—1)D’, | 
y =y—mD', 
2’ =z—nD', 


When these values of x’, y’ and 2’ are substituted in (30) and wx, J, m, n, are 
expressed in terms of y and z the result must be an identity in y and g, | 


weet ; 1 
A) —(l- 1)D = iu (y —mD yi + 25" (z —nD') + ONG = mD"')8 
+ Cisly—mD'\(z—nD'? + (a +C; ! 
& Pe + Cg ) (vy —mD') 
1 ‘ 
+ 4p," +Ci) (y—mD')?(z —nD' 


1 5 
In (ss +4) (z—nD')*+ vee Vee SS Oe Pee ee eres (33) 
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When the coefficients of the various powers of y and z in this identity are 
separately equated to zero, a set of relations between the primed and unprimed 
aberration coefficients is obtained, and these give the required transfer formulae. 

The direction cosines (/, m, n) are obtained as in § 6, but the values are required 
to a higher order of accuracy; equation (33) must be correct to O(p*), where p 
is of the same order of magnitude as y and x; hence / must be correct to O(p*) 
and m and n to O(p%), since D’ is O(p°). It is found that 


1 1 3 ine 4 9 Abea 
fet {aa T 73 bx + —- 7 C209” mG +2) Cupat (FC + 5 C30 + sa) 2" 


| 
eras | 1 1 | 
+ (2 (7 + ;) Coo +3Ca9Cio +2C?, =e 2B =a a) yrs? | 
we 4a 1¢C2 1 of O( p* : 
5 orate 2t 94 + O(p°); & 
| 
1 1 i 
a= {ry + 3C soy? + Cyo2” + 4C yoy? + (22 ae a) ya*h ee | 
| 


1 
vies F +2Cyoyz +2Cooy*s +4Cy2°| OD): 


Hence, since t=t' —D’, s=s’ —D’, it follows that 
y—mD' = ~Y+3D'Cyoy? + D'Cype? +4D' Coy? 
1 La 
i! oe aE yoy2 4 
+D (2¢e ae a) ya? + O(p4), 
z—nD' = ~ #+2D' Cz + 2D! Cyoy?2 +4D! Coys +O(p*). 


On substituting these and the value of x from (2) in (33) the identity in y and z 
is obtained; equating the coefficients of y? and 2? to zero gives the obvious 
relations 

t= =D"; eS e Oz) Oe Bree os (35) 


From the coefficients of y® and yz? the following are obtained, after some 
simplification : 


From the coefficients of y+, y?z? and +, after simplification including the use 


of (36) and (37), 


; BNE aie DDE. 5; 
Cy = (7) {Cio+ er ~ ci, Peo: (38) 
we LD gee Ds eas) 
Cos = (= = {Cia + a t Cie Cs5 + YZ Eee, Cy f > fF © ewe (39) 
s ; IB seco 
oe (=) {e% eee ay pi tee (40) 


Equations (36) to (40) are the transfer formulae for the aberration coefficients. 
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§8. DISCUSSION OF THE RESULTS 


The refraction and transfer formulae may be collected together as follows : 


Coefficient C30 Cie C40 C22 Cos 
Refraction formula equation no. (13) (28) (14) (29) (20) 
Transfer formula equation no. (36) (37) (38) (39) (40) 


It will be seen that the formulae for the fourth power coefficients Cyy, Cog, Cog 
contain one or both of the third power coefficients Cy), Cy. In the case of the 
refraction formulae this means that the amount of fourth power aberration intro- 
duced on refraction at a surface depends on the amount of third power aberration 
present in the incident wave front*. 

From the transfer formulae it can be seen that the wave front aberrations 
change in type as the wave front travels along the principal ray. For example, 
equation (40) shows that a wave front having primary coma (Cj,) will, after 
travelling some distance, also have some primary spherical aberration (Cg). 
This effect is usually quite small, but it may reach considerable proportions in a 
system with large separations between surfaces. By a similar argument to that 
in the last paragraph it can be shown that the change on transfer of an mth degree 
coefficient (C,;, where +7 =m) involves coefficients of the m— 1th degree. 

The formulae as given above are probably in the best form for exhibiting the 
relationships of the aberration coefficieats to each other and to the data of the 
principal ray trace, and in this paper they will be leftinthisform. Inasubsequent 
paper they will be put into forms more suitable for numerical computation and 
numerical examples will be given to show how closely the aberration coefficients 
reproduce the actual aberrations of a system, as determined by ray-tracing. 
Further work to be published includes formulae showing the effect of shifting 
the stop position, generalizations of the formulae for non-spherical surfaces, and 
formulae for the chromatic variations of the aberrations. 
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* There is a similar effect for all higher degree aberration terms. Suppose that the aberration 
coefficient C;; from the term y'z/ in (2) is being investigated, where i+j>4, so that quantities 
which are O(p' 4) must be retained in the expression for g ; since q is usually O(p), (unless i=0, 
when # is taken as 7/2 as in §5 and q is O(p2)) equation (5a) shows that terms in / which are 
O(0'*~1) must be retained, and from the method of evaluating J (see §5) it is clear that coefficients 
of terms in x which are O(p’*!—+) will enter into such terms in 1. That is, in the refraction 
formula for Cipp coefficients of degree one lower, such as Cea Cid j—g ete., will occur. In the 
case where 1=0, so that q is O(p?), this is still true, since y is O(p?) and the term in x containing 
C, fee in equation (5a) is O(p’); thus Cy jee occurs in the coefficient of p in the final expression 


for q. For example, if Cog were being evaluated, C,, would be nvolved, and so on. 
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The Temperature Variation of the Long-Wave Limit of 
Infra-Red Photoconductivity in Lead Sulphide and 
Similar Substances 


By=?t..S; MOSS* 


Physics Department, Telecommunications Research Establishment, Malvern 
Communicated by R. A. Smith; MS. received 13th May 1949 


-ABSTRACT. Measurements have been carried out on the spectral distribution of the 
photoconductive effect in layers of lead sulphide, selenide and telluride, at temperatures 
ranging from room temperature to that of liquid hydrogen. It is shown that cooling has 
a marked effect on the long-wave limit of sensitivity, which increases by as much as two 
“microns on cooling with liquid hydrogen. 
Brief descriptions of spectral measurements on these materials have already been published 
-by the author. It is shown that depopulation of energy levels due to thermal activation of 
-electrons is not responsible for the shift of spectral sensitivity. "The suggestion is made that 
the shift may be related to the variation of dielectric constant with temperature. 


§1. INTRODUCTION 
CONSIDERABLE number of materials, mainly oxides and sulphides or 
A related types, are now known which exhibit photoconductivity on 
irradiation by light in the visible or near infra-red region. Examples 
are HgS, CdS, CdSe, CdTe, S, Se, MoS,, Sb.S,, Bi,S,, ZnS, Ag,S, Cu,O, 
PbS, PbSe, PbTe, PbCrO,, and a few mixed crystals. 

For stibnite, Sb,Ss3, Elliot found that the wavelength of peak sensitivity was 
reduced from 0:76 at —20°c. to 0-68 at —190°c. (Elliot 1915). 

For acanthite, Ag,S (Coblentz 1919), it was found that the wavelength of 
peak sensitivity decreased from 1-5 to 1-34 on cooling from 22°c. to —158°c. 
Similar results were obtained for argentite—an alternative crystalline form 
of Ag,S. e%. 

For stephanite (Barish 1935) a progressive reduction in the threshold wave- 
length from 0-8 » to 0-74 takes place as the temperature is reduced from 21° c. 
to —125°c. par 

The absorption edge of MoS, moves from 0-7 x, at room temperature to 0-66 wu 
at —190°c. (Coblentz 1919). 

Frerichs (1947) showed that for CdS there was a linear relation between 
the temperature and the threshold wavelength which varied from 0-58 uw at 500° c. 
to 0-49 w at —183°c. 

For the substances PbS, PbSe, PbTe, which show conductivity at longer 
wavelengths than the substances mentioned above, the threshold wavelength 
increases with cooling. Experiments on this effect are described in the following 
section. 

§2. EXPERIMENTAL RESULTS 

The material used for the measurements of spectral sensitivity was in the 
form of microcrystalline layers deposited on the inner walls of small evacuated 
Dewar vessels. The PbS layers used were chemically deposited, and sub- 
-sequently sensitized by heating in oxygen. 


* Now attached to the Department of Colloid Science, Cambridge. 
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The PbSe and PbTe layers were made by evaporation in vacuo of appropriate 
quantities of the bulk compounds. ‘These compounds were prepared by heating 
the constituent elements together in evacuated quartz tubes for some hours 
at 1000°c. Both the PbSe and PbTe were made with an initial excess of Se 
(or Te) of 5% over the stoichiometric composition, and then purified by vacuum 
sublimation. 

The material was sensitized by heating at low pressure (about 10-°mm. Hg) 
at a temperature just below the sublimation temperature for the compound, for 
periods varying from one to twenty hours. The final layer was then formed by 
evaporating the film by a blowtorch flame whilst the electrode surface was water- 
cooled. 

It was necessary to provide the cells with windows capable of transmitting 
the range of wavelengths required. For PbS, the window was an “‘infra-red glass” 
of the soda—zinc-oxide-silicate type. The transmission of the window was 
measured after the other measurements on the layer were completed. It was 
found that the transmission was more than 70°, between 1-0 and 4-2, falling 
to 25% at 45y. The spectral sensitivity curves were corrected for this loss. 

For the PbSe and PbTe cells, the window material was KRS 5 (a mixture of 
thallium bromide and iodide) which had 65° transmission from 1 to 35. Plates 
of this material were fixed to the cell envelopes by vacuum wax. Thus the 
envelopes could not be heated during the sensitizing process. "Two means of 
overcoming this difficulty were evolved: 

(a) The window was kept cool by immersion in water, whilst the layer 
(deposited on the inner wall) was heated by filling the inside of the cell with a 


low melting point alloy, maintained at the requisite temperature of about 400° c. 


by a small electric heater. 

(5) ‘The envelope was made with an elongated neck to which the window was 
attached. ‘The envelope was then mounted on the vacuum system with the window 
protruding through the base of an oven, so that it could again be watercooled. 

No measurements of the reflection or absorption of the layers were made. 
However, it is known (Mellor 1923) that the reflection coefficient of PbS is constant 
at 31%, over the wavelength range of 1 to 10u. For PbSe the reflection coeffi- 
cient for red light is 40°% (Schneiderhohn 1931) and it will probably stay constant 
at about this value in a like manner to that of PbS. Similar behaviour would 
be expected for PbTe. ; 

To overcome the difficulty of the absorption varying with wavelength, the 
layers were made relatively thick, so that there would be almost no transmission 
at all wavelengths used. ‘Then, with the reflection coefficient roughly constant, 
the energy absorbed would be proportional to the incident radiation. 

Measurements on PbS films* at room temperature have shown that the 
absorption varies only gradually from 1p to 4. Similar results have been 
obtained (Simpson and Daly, unpublished data) for PbTe, where a layer about 
0-2 thick had 80°% absorption in the near infra-red, falling to 50% at longer 
wavelengths. 

The thicknesses of the layers used by the author were estimated by weighing 
the films after the rest of the measurements were completed. The values 
obtained were: PbS 1-8y, PbSe 1-24, PbTe 1-0 (approximately). 

The spectrometer used in these measurements was a double monochromator 
with lithium fluoride prisms and having a high degree of freedom from scatter. 

* Unpublished experiments carried out by E. J. Harris at T.R.E. in 1945. 
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Radiation from a Nernst source was interrupted by a motor driven chopper disc 
at 800 c/s. before entering the spectrometer. The resultant signal from the 
cell was fed into a high gain amplifier tuned to 800 c/s., the output from which 
was read on an A.c. meter. 

A Hilger thermopile was used to measure the energy emerging from the- 
spectrometer. 

The sensitivity figures for an equal-energy spectrum could then readily be 
found by virtue of the fact that, for the intensities used, the signal from the layer: 


i 
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Figure 1. Spectral sensitivity of lead sulphide. Scales adjusted to make peak sensitivities equal. 


was proportional to the incident radiation. It may be noted that, due to the 
manufacturing difficulties caused by the use of windows which could not be 
heated, these layers were not of optimum sensitivity. For such layers the 
linearity law holds up to higher intensities. 


§3. MEASUREMENTS ON PbS, PbSe, PbTe 
For lead sulphide, measurements were carried out at room temperature, 
and when cooled by solid carbon dioxide in acetone, liquid oxygen and liquid 
hydrogen. ‘The results (referred to an equal-energy spectrum) are shown in 
Figure 1. 
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It will be observed that the curves are comparatively flat at the shorter wave- 
lengths, but that as the wavelength is increased, a “ threshold Is reached beyond 
which the sensitivity decays rapidly in an approximately exponential manner. 
“This point extends to longer wavelengths as the temperature is lowered. Defining 
the “threshold” as the wavelength at which the sensitivity has fallen to 50% of 
its peak value, we see that cooling has increased this wavelength from 2-9 to 
4-2 microns. ; s 

For the PbSe layers sensitivity measurements were made at 290) Key ho) ee 
and 90°x. The room temperature measurements were made on a different layer, 
since the one used for the 90°xK. and 195°x. results had insufficient sensitivity 
when uncooled. The results obtained, which are shown in Figure 2, resemble 
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Figure 2. Spectral sensitivity of lead selenide. 


those for the sulphide. As the curves are flat over a considerable range of wave- 
lengths, it is convenient to take the threshold where the sensitivity has dropped 
to 50% of its value at some point on the flat part of the curve, say 3-04. With 
this: definition, it is seen that the threshold is extended from about 3-5 to 5-0u 
on cooling to 90° K. 

For PbTe, experiments were carried out at temperatures of 20°K., 77°K., 
90°K., 195°K.and234°xK. Thislatter temperature, the melting point of mercury, is 
the highest temperature at which any of the PbTe layers yet made had adequate 
sensitivity to enable a spectral measurement to be made. Even so, the sensitivity 
at both 234°K. and 195°k. was poor, with consequent loss of accuracy—particu- 
larly at the longer wavelengths. ‘The experimental results are plotted in Figure 3. 
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The measurements at 20°. did not cover the complete wavelength range- 
due to the fact that the design of the cell envelope was unsatisfactory for use with 
liquid hydrogen, which evaporated so rapidly that only a small number of points. 
could be plotted. These were accordingly taken over the threshold region. 
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Figure 3. Spectral sensitivity of lead telluride. 


As before the curves are very flat, and the threshold wavelength is again taken: 
where the sensitivity is half that over the flat region. It will be seen that this 
wavelength is increased by cooling from below 4:0 to 6-0 u at 20°K. 


§4. DISCUSSION OF RESULTS 


It has been found for all three types of layer, that if the quantum value 
appropriate to the threshold wavelength is plotted against the temperature, the 
points fall on a straight line. Such plots are shown in Figure 4. The PbS 
results are shown on the top line, while those for PbSe and PbTe lie on the 
same (lower) line. 

The threshold quantum value «, may be expressed as ep=€)+a%7, where 
a has a similar value for all three materials. Any of the family of PbS spectral 
sensitivity curves can be represented by an empirical formula of the type 
S=A/[1+Be”“] where 2 is in microns. a is determined from the slope 
of the curve beyond the “‘threshold’’, B is adjusted to give the threshold at the: 
correct wavelength, and A is merely a scale factor for the absolute sensitivity. 
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B is independent of temperature and has the value e, while a@ varies with 
‘temperature as follows: 
TECK.) 290 Mee) 90 20 : 
a 9:92 8-90 7:62 6:85 
From the above relation it may be seen that S will be half its short-wavelength 


value when 


Be“ =1, 
i.e. the threshold wavelength Ay is given by 
the | 


Since B is constant, then ax1/Ap«e, giving a linear relation between a and 
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Figure 4. 


‘the temperature which is found to be a=(575+T)/87. Thus the whole family 
of spectral sensitivity curves for PbS may be expressed as 
Ha SUES, 
1 + exp [A(575 + T)/87 —29] 
‘This function is plotted in Figure 5 for comparison with the experimental curves. 
The PbSe and PbTe results may be expressed by a similar formula, i.e. 
epee i eee al 
1+ exp [A(388 + T)/125 — 19] 


In this case, however, the agreement is inferior to that for PbS. 


§5. POSSIBLE CAUSES OF SPECTRAL SHIFT 
So far there is no general theory of photoconductivity in semiconductors, and 
the factors determining the threshold wavelength and its temperature dependent: 
are not yet fully understood. 
One factor which it was thought might contribute to the shift of threshold 
wavelength is the varying density of conduction electrons in the material as the 
temperature changes. On the band theory of solids, increasing numbers of 
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electrons will be raised from the full band to the conduction band of energy levels 
as the temperature is raised. It is possible that these electrons could fill the lower 
levels in the conduction band to such an extent that the optical transitions giving 
rise to the photoconductivity could only take place to levels somewhat above the 
bottom of the conduction band. Depopulation of the levels at the top of the filled 
band would similarly modify the optical transitions. The effect would be 
emphasized if the edges of the energy bands were not sharply defined, but tailed 
off slowly, so that a relatively small number of levels would have to be filled to. 
produce the change of 0-1 ev. required to shift the threshold wavelength the 
observed amount. Soft x-ray work (Skinner 1938) indicates that there are 
‘tails to the bands. 
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Figure 5. 


The mechanism suggested could be checked experimentally by using a strong 
light to produce the high density of conduction electrons, whilst maintaining the 
temperature constant. Accordingly, the spectral sensitivity was measured for 
a PbS layer at room temperature and then at liquid oxygen temperature, with a 
steady illumination on the layer such that the conductivity was almost restored 
to its original room temperature value. It was found that the usual large shift 
occurred on cooling, and that the application of the light had a negligible effect. 

Another mechanism which may contribute to the shift of the spectral sensi- 
tivity curves is change of the refractive index of the material. If the photo-effect 
does not take place equally at all lattice points but occurs only at privileged 
positions where the periodicity of the lattice is disturbed by either interstitial 
ions, vacant lattice sites, or cracks, then these centres will be somewhat similar to 
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isolated atoms except that they are immersed in a medium of dielectric constant 
equal to that of the bulk material. The energy levels of the centres will then be 


reduced in the ratio 1/K®, where K is the effective dielectric constant of the 
material in the neighbourhood of the centre. The threshold of the photo- 
conductive effect will then depend on the dielectric constant, any temperature 
variation of which will cause a spectral shift. 

According to Mott and Gurney (1940) this dielectric constant will be nearly 
that of the material at optical frequencies—i.e. (refractive index)’. 

Unfortunately no data are at present available on the temperature variation 
of the refractive index of the PbS series with which to test this hypothesis. * 


§6. CONCLUSIONS 


Results are given of measurements carried out on the spectral distribution 
of infra-red photoconductivity in lead sulphide, selenide, and telluride layers at 
various temperatures. It is found that the long-wave threshold of sensitivity 
moves to longer wavelengths as the temperature is reduced, the quantum value 
of the threshold being a linear function of the absolute temperature. No full 
explanation of these results is yet available, but it is suggested that temperature 
variation of the refractive indices of the materials may be the cause of the shift 
of spectral sensitivity with temperature. 
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REVIEWS OF BOOKS 


Mercury Arcs, by F. J. Teaco and J. F. Git. Pp. viii+107. Methuen’s 
Monographs on Physical Subjects. 2nd (revised) Edition. (London: 
Methuen, 1949.) 6s. 6d. 


The mercury arc rectifier is the most common form of gas-filled valve, and is of great 
interest to the modern electrical engineer. The control of very high currents, voltage 
regulation, frequency changers and speed control of heavy a.c. motors, are just a few of 
the modern applications which have now developed into an important branch of electrical 
engineering. 

This type of rectifier operates in virtue of the copious production of electrons at local 
regions of the cathode, which is a pool of mercury; the electrons produced ionize the vapour 
through which a low-voltage arc passes to the anode. Voltage drop is practically inde- 
pendent of the current passed. Vapour is supplied by the evaporation of the cathode, and 
mercury has the advantage that this loss due to evaporation is made good by the return of the 
liquid condensed from the vapour at the cool walls of the containing vessel. If no electrons 
can. be produced at the other electrode (the anode) then no current passes when the applied 
voltage is reversed, and perfect rectification takes place. On the other hand, if for some 
reason the anode develops electron-emitting regions, then a glow discharge, which possibly 
leads to an arc, can take place on the reverse phase, and rectification is reduced ; when a full 
arc (back-fire) occurs there is complete failure of the rectifying action. Thus the mercury 
arc rectifier depends for its action on the asymmetry of the properties of the cathode and 
anode. 

Now the operation of such a rectifier presents a number of fundamental problems of 
great interest to the physicist. In the first place, there is the important problem of the 
mechanism of electron production at the cathode hot-spot, and various mechanisms have 
been suggested :—local thermionic emission ; cold emission in the field set up by the positive 
ion space charge near the cathde hot-spot ; microscopic oscillation of the mercury surface, 
set up by fluctuations in the surface, leading to sharp projections of mercury from which 
emission then occurs due to the intensified field ; but the full explanation has yet to be given. 
In fact, considerable doubt exists concerning the actual state of matter in the microscopic 
region near the surface of the hot-spot. Moreover, the development of electron emitting 
regions in the anode presents an interesting problem to the physicist. Consequently, a 
modern book on mercury arcs is likely to be read with great interest by the research physicist 
and electrical engineer alike, to see what discussion is given to these fundamental problems. 

Mercury Arcs by Professor F. J. Teago and the late J. F. Gill was first published in 1936, 
but a second edition, revised, has been issued this year. It has the object of outlining the 
“more important aspects of the arc from the electrical, mechanical, and operation stand- 
points ; ....and fine detail has purposely been omitted in order that the essential principles 
may stand out clearly’. 

The book has 11 chapters and an appendix. The first chapter, of 12 pages, covers the 
theory of the rectifier, including cathode processes and back-fires. Chapter 2, of 15 pages, 
gives constructional details of a high-output tank rectifier, together with some observations on 
high vacua, including diffusion pumps, McLeod and Pirani gauges. The rest of this book 
of 104 pages is taken up with what may be called circuit analysis, as indicated by the chapter 
titles : Grid control of voltage ; Grid control of current ; Voltage regulation ; Current- 
voltage characteristics ; Rectification, inversion and regeneration ; Wave-form analysis ; 
Transformer ratings ; Power factor ; Mathematics of mercury arcs ; Appendix, The 
Pirani gauge. The chapter on the mathematics of mercury arcs has nothing at all to do with 
the mechanism of the arc itself ; it does not discuss the diffusion of ions and electrons, 
recombination, and ion-bombardment of electrodes, but gives analyses of transformer 
circuits containing ideal rectifiers. 

It should thus be clear that the title of the book is misleading ; a more appropriate title 
would have been Mercury Arc Circuit Analysis. Considered as a book on Mercury Arcs, 
it fails in the aim stated in the Preface, viz. to make essential principles stand out clearly. 
The inadequate treatment of the physics of mercury arcs is the more unfortunate as the 
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book appears in Methuen’s Series of Monographs on Physical Subjects. The brief dis- 
cussion given on the mechanism of the arc, which is the subject of greatest interest to the 
physicist, contains some questionable statements. For instance, the rectifying action does 
not as is stated (page 1) depend upon the unilateral conductivity of an ionized vapour. 
In the discussion (pages 4 to 6) of the cathode hot-spot, the impression is given that electron 
production at the cathode is accounted for by thermionic emission, although experiments 
are quoted from which the authors conclude (page 8) that it is doubtful if the hot-spot 
temperature ever reaches 1000°c. Again, the mechanism of back-fire is treated very 
briefly for so important a subject, and the practical problem of commutation oscillations is 
not dealt with at all. Considering the inadequate treatment of these fundamental problems 
of the mercury arc rectifier, it is difficult to find justification for the inclusion of an appendix 
entitled the Pirani Gauge, which in fact gives the well-known theory of the unbalanced 
Wheatstone bridge and nothing else. It appears that the book is likely to be of more interest 
to the circuit engineer than to the physicist. F, LLEWELLYN JONES. 
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ABSTRACTS FOR SECTION A 


On the Size-Frequency Distribution of Penetrating Showers, by W. HerrLer and 
L. JANossy. . 


ABSTRACT. We have calculated the relative number of showers containing » mesons 
when a fast primary nucleon hits a compound nucleus on the assumption that in a single 
nucleon-nucleon collision only one meson is produced (pure plural production). The 
cross section for meson production is assumed to have the form ®(¢/E) de/E, where E is the 
energy of the primary and « the energy lost by it, but the absolute value of ® and its depen- 
dence on ¢/E are left open. The fluctuations of the shower size are duly taken into account. 
For comparison with experiments it is assumed that about half the number of fast shower 
particles are fast protons as suggested by the positive excess of penetrating shower particles. 
The theoretical results depend on two constants, the total cross section and some weighted 
average over the energy distribution O(</E). Choosing these suitably, the results are in 
good agreement with the relevant data within the experimental statistical errors, and the 
agreement extends also to the “ large showers ”’ with 20 or so thin tracks. The values of 
the constants thus determined are physically reasonable. It is concluded that the size— 
frequency distribution of penetrating showers can be fully understood on the basis of pure 
plural production, but some—small or large—contribution from multiple processes cannot 


be excluded. 


The Meson Intensity at the Surface of the Earth and the Temperature at the 
Production Level, by A. DUPERIER. 


! ABSTRACT. Meson intensity has been recorded by using a triple-coincidence arrangement 
) with a lead absorber 25 cm. thick between the counters. The comparison of daily mean 
i values with the height of a number of pressure-levels has led us to take into consideration 
| the temperature of the air layer between 200 and 100 mb. The results show that this 
% temperature appears to be a factor closely controlling the meson intensity at the surface 
of the earth; the intensity increases with increasing temperature at the rate of 0-:12% per °c. 
Other factors in determining the intensity are the mass of air and the height of the 100 mb. 
level over the station. The corresponding coefficients of absorption and decay prove to be 
>» comparable with those obtained from other measurements for mesons of the same momentum. 
‘y It has also been found that the coexistence of the three effects permits explanation of the 
variability exhibited by the barometric coefficient when evaluated from short periods of 
observations. 
These results may therefore be taken as evidence that the bulk of mesons is produced 
7 at the height of 100 mb. 
| If the positive effect of temperature is interpreted as the result of processes of decay and - 
| interactions with air nuclei of 7-mesons at the production level, then the coefficient 
0-12% per °c. leads to a value of 4:9 x107® sec. for the upper limit of the mean life of 
m-mesons. 


On the Decay of 7-Mesons, by C. B. VAN WYK. ~ 


1 ABSTRACT. The production of t-mesons (mass about 900 electron masses) by cosmic 
} rays in rather thin layers of air in the upper atmosphere seems to indicate that the nucleon 
} T-meson coupling constant is not very small. ‘The Berkeley experiments show that this is 
4) also true for the nucleon—7-meson coupling constant. Consequently the decay of a charged 
' ¢-meson into a photon and a charged 7-meson by means of virtual creation and annihilation 
of nucleons ought to be fairly probable. In this paper this decay process, with the emission 
of one photon, is studied for the various types of r- andz-mesons. Unless both mesons have 
_ zero spin, the matrix elements, and therefore the lifetimes, are found to be of the same order 
of magnitude, and the r-meson turns out to be too short-lived to be observed photo- 
4 graphically. 
If both mesons have zero spin the emission of only one photon is forbidden. The 
| detailed discussion of higher order decay processes in this case is reserved for a later paper. 
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Excited Electronic Levels in Conjugated Molecules—III: Energy States of 
Naphthalene, by J. JACOBS 


ABSTRACT. ‘The electronic states of naphthalene are determined using antisymmetrized 
molecular orbitals, which are based on wave functions given by the simple molecular orbital 
theory. This method allows the electronic interaction term to be included explicitly in the 
Hamiltonian for the system. A clear distinction is made between electronic configurations 
determined with these wave functions and the actual states of the molecule: states arise 
from a superposition of configurations of the same symmetry. The ground state becomes 
considerably stabilized by this configurational interaction, to which even highly excited 
levels contribute. Calculated energy values for some of the lower tiansitions agree 
well with ultra-violet spectral data, as do the newly derived oscillator strengths. A tentative 
assignment is made regarding the polarization of the electronic transitions in the region of 


30,000—55,000 cm™?. 


A New Technique for the Spectroscopic Examination of Flames at Normal 
Pressures, by H. G. WoLFHARD and W. G. PARKER. 


ABSTRACT. Up to the present spectroscopic studies of flames supported on a tube 
have been limited by the geometry of the flame. Thus the important reaction zone presents 
a three-dimensional problem in which it is impossible to locate the emitters exactly. 
Moreover this zone is so thin in premixed gas flames that intermediate products cannot be 
determined from absorption studies. "To overcome these difficulties the authors have 
constructed a burner with a two-dimensional diffusion flame in which the reaction zone at 
ordinary pressures is 5-10 mm. thick and has an optical depth of 5 cm. or more. With 
this flame many of the reacting molecules can be located through either their emission or 
absorption spectra or both. The value of the new technique in the elucidation of com- 


bustion processes is demonstrated with reference to NH;—O, and H,—Oy flames, and mention ~ 


is made of its application to hydrocarbon flames. 


Dielectric Changes in Phosphors containing more than One Activator, by G. F. J. 
Garuick and A. F. Gipson. 


ABSTRACT. It is well known that some luminescent materials, notably zinc sulphide, 
show marked changes in dielectric constant and dielectric loss when excited. Recent work 
by the authors has established that the changes are due to the filling of electron traps, which 
have large effective diameter and high polarizability. 

The present paper describes measurements of the dielectric changes in phosphors 
containing two or more luminescence activators. Assuming that the dielectric changes are 
due to electron trapping, it is possible to derive information regarding the nature of electron 
traps and in particular their apparent association with the luminescence centres. The 
results support the contention that the introduction of the activating impurity into phosphors 
produces trapping states for excited electrons and that each trap forms part of a complex 
which includes the luminescence centre. 


Eigenfunctions following from Sommerfeld’s Polynomial Method, by A. RUBINOWICZ. 


ABSTRACT. The eigenfunctions following from Sommerfeld’s polynomial method were 
determined in the general case. They are products of a convergence factor E and a poly- 
nomial P given by a hypergeometric series. 'The form of E depends upon whether P is 
expressible by the ordinary or the confluent hypergeometric series. 


CORRIGENDUM 
“ The Physical Basis of Life”, by J. D. Bernat (Proc. Phys. Soc. B, 1949, 62; 597). 


Page 600, last clause of footnote , viz.* : ‘‘ and secondly, they must originate and develop 
out of some pre-existing system, or, in plain English, they must work and they must have 
got there in the first place ” should occur in the text immediately after the words “ dynamic 
stability ’’ (line 22). 
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